


FEASIBILITY STUDY FOR A 
MICROWAVE-POWERED 
OZONE SNIFFER AIRCRAFT 


VOLUME I 


M r L : ! i'/L-i 
. . T i ; ' i :> 

J h 


'-u) t 1 A'- I H ITY Ji T u '' Y 
■ y ? ... L i * ' A ' f r K A T hi 
f vi n r r > t * r i > 1 v l ^ c ■ ’ n i C 


1 Hk a 


. : /"> , , 
j'ti ■’ y , j 

K A i- T 



ini>U ' 



C L * ’ i V - 


Unc 1 


• /u r . 

u 


; ' / 


FEASIBILITY STUDY FOR A 
MICROWAVE— LOWERED 

OZONE SNIFFER AIRCRAFT 

Project #: MSG-89A1 

Major Qualifying Project Report 
Submitted to the Faculty 
of the 

Worcester Polytechnic Institute 

in partial fulfillment of the requirements 

for the 

Degree of Bachelor of Science 


by 

David F. Botros 





Charlotte K. Co 


w 


Noah P. Forden 



Martin A. Hels^ffg 

j2a IL 

Thomas H. Jutras 


DoHoon Kim 



_ cPjL 


Ethan M. Odin 

\>- - 

Scott B. Sandler 


Approved By: 


Melvin S. Garelick, Asst. Professor 


Joseph J. Rencis, Professor 



Abstract 


The preliminary design of a high-altitude, remotely- 
piloted, atmospheric-sampling aircraft powered by 
microwave energy beamed from ground-based antenna has 
been completed. The vehicle has a gross weight of 6720 
pounds and is sized to carry a 1000 pound payload at an 
altitude of 100,000 feet. The underside of the wing 
serves as the surface of a rectenna designed to receive 
microwave energy at a power density of 700 watts per 
square meter and the wing has a planform area of 3634 
square feet to absorb the required power at an optimum 
Mach number M = 0.44. The aircraft utilizes a horizontal 
tail and a canard for longitudinal control and to enhance 
the structural rigidity of the twin fuselage 
configuration. The wing structure is designed to 
withstand a gust-induced load factor n = 3 at cruise 
altitude but the low-wing loading of the aircraft makes 
it very sensitive to gusts at low altitudes, which may 
induce load factors in excess of 20. A structural load 
alleviation system is therefore proposed to limit actual 
loads to the designed structural limit. Losses will 
require transmitted power on the order of megawatts to be 
radiated to the aircraft from the ground station, 
presenting environmental problems. Since the 
transmitting antenna would have a diameter of several 
hundred feet, it would not be readily transportable, so 
we propose that a single antenna be constructed at a site 
from which the aircraft is flown. The aircraft would be 
towed aloft to an initial altitude at which the microwave 
power would be utilized. The aircraft would climb to 
cruise altitude in a spiral flight path and orbit the 
transmitter in a gentle turn. 
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1.0 introduction 


NASA requires an unmanned flight vehicle which can 
obtain atmospheric samples at an altitude of 100,000 
feet to determine the level of constituents such as 
ozone . 

Recent reports about the rapid depletion of the 
ozone layer have caused much concern among 
environmentalists worldwide. The ozone layer protects 
the earth from the harmful radiation of the sun. 
Preliminary information indicates that fluorocarbons 
emitted by aerosol products and refrigeration equipment 
have caused a hole in the ozone layer over the south 
pole. More recent studies also show signs of ozone 

depletion around the equator. 

Currently there is no method of obtaining accurate 
atmospheric samples of ozone. Conventional aircraft are 
unable to operate at the extreme altitudes required for 
ozone sampling, due the lack of air required for fuel 
oxidation in an air-breathing engine. Also, the low air 
- density at altitude causes difficulty in providing 
sufficient lift. Rockets are able to reach these 
altitudes, but only for a short duration and at high 
velocities. In order to make accurate measurements of 
ozone concentrations the atmosphere must be undisturbed 
by the aircraft. 

The unique requirements of this mission necessitate 
the design of a long-endurance, high-altitude aircraft. 
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2.0 Mission Profile 


An aircraft is required to fly an atmospheric 
sampling payload weighing 1,000 lbs. to an altitude of 
100,000 feet and loiter for extended periods, possibly 
on the order of weeks. The high altitude and extended 
duration of this mission makes a manned vehicle 
impractical, as life support systems would add a great 
deal of weight and complexity. Therefore, a remotely 
piloted design is necessary for the completion of the 
mission . 

The aircraft is equipped with a microwave power 
system, which receives energy beamed to the aircraft 
from a single transmitting station on the ground. This 
power can be transmitted for indefinite periods of time, 
and allows the vehicle to remain airborne for as long as 
the mission requires. The microwave power is 
transmitted in a conical pattern, so the flight path is 
constrained to a spiralling climb to altitude. A 
diagram of the flight pattern is shown in Figure 
2 . 1 . 

The aircraft is designed to be towed to an altitude 
of 20,000 feet, at which point the entire flight path is 
contained within the conical microwave beam. With a 
maximum climb rate of 11.88 ft/sec, it attains cruise 
altitude within four hours of takeoff. Once at the 
design altitude of 100,000 feet, the flight path radius 
is approximately 28,500 feet. The aircraft has the 
capability to loiter in this condition almost 
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indefinitely. The cruise speed is at a Mach number of 
0.44, which allows operation well below velocities which 
would induce areas of transonic flow over the aircraft. 
Atmospheric sampling may then be accomplished before 
descent. 
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3.0 Propulsion 


3.1 Power Systems 


Various types of propulsion system that could be 
used to power the aircraft were considered. The design 
requirement for a vehicle capable of flight at 100,000 
feet for an extended period of time limited the choice 
of propulsion systems to solar, hydrazine, diesel, 
propane, regenerative fuel cells and microwave. A weight 
estimation program was used to determine the range of 
power required (see chapter 4). 

Light weight diesels will produce peak power at a 
maximum altitude of only 70,000 feet, which does not 
meet mission equirements. 

Another possible power source was hydrazine. One 
aircraft design used hydrazine to reach a projected 
altitude of 100,000 feet (see Ref. 7) but the endurance 
is limited. The same result was observed with propane 
powered flight. 

One attractive aspect of solar power is the large 
planform of the wing that can be used for the location 
of the photo-voltaic cells. Using thin film technology, 
the vehicle would not have the same altitude limitations 
as air breathing engines and would also benefit from a 
low power to weight ratio. Solar power has a projected 
conversion efficiency of approximately 15%. which 
results in an expected power density of about 195 watts 
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per square meter. When the weight estimation program was 
run using a wing area proportional to this power 
density, the program would not converge on a solution. 

It was concluded that solar energy would not meet the 
initial requirements. Solar power would only provide 
about 15% of the total power requirement of the vehicle 

during peak hours of sunlight. 

The concept of using microwave energy to power an 
aircraft involves producing the power at a distance away 
from the aircraft and then beaming the energy in the 
form of microwave radiation to the aircraft. The vehicle 
has a receiving antenna which has rectifying circuitry 
that converts microwave power into DC power. The current 
produced is used in an electric motor that drives the 
propellers. The nature of the system is such that the 
power obtained is directly proportional to the size of 
the rectenna. A power density of 700 watts per square 
meter is presently considered attainable. The microwave 
energy system weighs approximately 0.79 kilograms per 
kilowatt of power needed. To maximize received power and 
minimize aerodynamic drag due to the rectenna, it was 
decided to place the rectenna on the underside of the 
wing (see Figure 3.1). This produced a direct 
relationship between the power density and wing area 
similar to the solar powered configuration. When the 
weight estimation program was run using the revised 
power density, the solution converged. 
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Another option was solar energy as an addition to 
the microwave power with both routed into a regenerative 
fuel cell. This would give additional power from the sun 
while it was available and the fuel cell would ensure 
continuous power supply, but this advantage did not 
compensate for the additional weight of the fuel cell. 
The most efficient cells available were hydrogen-oxygen 
regenerative fuel cells, however, the system is only 
capable of 1\ horsepower, while weighing 500 pounds. 
Since this power represents only 2 to 5% of the required 
power for the plane, it became clear that this approach 
was not feasible. 

3.2 Microwave Propulsion 

With the elimination of all other alternatives, 
microwave energy was preferred. With the power source 
decided, the drive train could be designed. 

Energy will be routed directly from the rectenna to 
the electric motor. The brushless, rare earth magnet 
electric motors used have a power density of 2 to 4 
horsepower per pound and an efficiency of approximately 
97%. A frequency of 2.45 gigahertz was selected using 
data on general power efficiency and degradation as a 
factor of climate and frequency (see Figure 3.2). 

were four possible methods of implementing 
the microwave propulsion system. First, it would be 
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possible to place a satellite or series of satellites in 
orbit and use massive solar panels to collect energy 
from the sun and beam it down onto the plane in flight. 
This idea was not pursued because satellite usage was 
deemed uneconomical. Another choice involved a second 
aircraft orbiting at a lower altitude with a large 
generator and transmitter onboard (see Figure 3.3). This 
plane would be kept on station using air-to-air 
refueling and rotate its schedule with other aircraft to 
provide an around the clock power source to the ozone- 
sniffer. This system of support aircraft would allow for 
redeployment of the research aircraft to any point in 
the world. This alternative was not selected since the 
size of the transmitting antenna (30 to 50 meters m 
diameter) would be difficult to accomodate. The third 
and fourth systems made use of semi-permanent ground 
stations. One consists of a chain of ground stations 
linked to provide a corridor in which the plane would 
fly (see Figure 3.4). The other used only a single 
source, forcing the vehicle to fly in a circular path to 

stay within the beam (see Figure 3.5). 

The last alternative of providing the microwave 
energy to the aircraft was selected because the ground 
stations will be large (approximately 100 meters in 
diameter) and expensive to build. When the aircraft is 
in use, it may not be possible to locate the ground 
stations at every location desired. With only one 
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station in use, the problem is not as severe. The steady 
turning flight, however, requires more power than level 
flight and also reduces the angle of incidence of the 
microwave beam on the rectenna, thereby reducing 
available power. This effect is examined in chapter 9 
and will not prevent the aircraft from flying. 

3 . 3 S umm ary of the Propulsion System 

The high altitude research aircraft will use 
microwave energy beamed up from a single ground station 
that it will circle above. The vehicle will convert the 
energy into DC power and use that to run an electric 
motor that will drive a propeller. Other types of 
propulsion system could be used, but the system chosen 
was selected as the best possible. 
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4.0 Aircraft Sizing 


In order to provide an initial component weight 
estimate, a computer code, WEIGHT. BAS, was developed. 
This program, written in BASIC for IBM compatible PC's, 
is based on the weight equations given in Reference 24 
for light utility aircraft. The sizing is done at the 
design altitude of 100,000 feet. The fact that the 
aircraft is microwave powered necessitated modification 
of the basic equations. The weight of the fuel, which 
is usually a significant portion of a conventional 
aircraft's gross takeoff weight, was omitted. In place 
of the fuel equation, an equation for the weight of the 
rectenna was derived from data in Reference 6. An 
additional modification in the form of a correction 
factor for composite materials has been applied to the 
final component weights of the wing, fuselage, 
horizontal and vertical tails. This correction factor 
is also based on data given in Reference 24. Once the 
component weights were calculated, an aircraft 
configuration was determined. 

4.1 Component Weight Estimation 

Several factors complicated the estimation process. 
First, there is little available data about sizing high 
altitude platforms on which to base an initial design. 
Most initial aircraft sizing is based on a comparison of 
mission requirements with the mission requirements of 
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existing aircraft. Without the advantage of hindsight, 
the initial sizing began by examining the variation of 
the wing loading (gross aircraft weight divided by the 
wing planform area) with the Mach number. This 
calculation is carried out over a range of lift 
coefficients from .1 to 1.2 as shown in Figure 4.1. At 
each Mach number the maximum subsonic lift coefficient 
was determined using the Kdrmdn— Tsien relation given on 
p. 258 of Reference 1. This subsonic lift coefficient 
is the maximum value obtainable without local velocities 
reaching Mach 1. 

Only flight velocities less than Mach 1 were 
considered for this mission, due to the increased power 
which would be required at transonic or supersonic 
speeds. Data from several studies of high altitude 
platform design (References 2,7,16,22) were plotted onto 
Figure 4.1. The design of Reference 7 was a high 
altitude, hydrazine powered aircraft while Reference 16 
was a microwave powered configuration. Particular note 
should be made that most designs with liquid fuel were 
located on the graph in the region of the Reference 7 
design. Most aircraft that did not have onboard fuel 
were in the range near the Reference 16 configuration. 

Since all of the component estimation equations in 
Reference 24 are based on statistical data for 
conventional aircraft, these equations may not 
accurately represent the microwave-powered 


14 



configuration. At low wing loadings, the wing weight 
equation is very sensitive to small changes in wing 
area. The change in wing weight with wing area, dW/dS, 
varies as approximately .61S~”. However, for lack of 
a better methodology and statistical equations, the 
procedure of Reference 24 has been used. 

For a microwave powered aircraft, the size of the 
wing limits the maximum power that the aircraft is 
capable of absorbing since the wing acts as the 
rectenna. However, a separate rectenna could be used 
to absorb any required power in excess of that which is 
received by the wing. The external rectenna would look 
similar to the radar dish on an AWACS aircraft except it 
would be mounted under the fuselage. 

4.2 Discussion of Parametric Study 

The series of graphs generated from the weight 
estimation program demonstrated many of the tradeoffs to 
be considered in the selection of the design parameters 
of the aircraft. These graphs are shown in Figures 4.2 - 
4.7. The output consisted of six graphs. Within each 
series are graphs of gross weight, wing span, wing area, 
lift-to-drag ratio, and external rectenna area vs. Mach 
number for lift coefficients of .5, .6, .7, .8, and .9. 

Superimposed on these graphs are lines denoting constant 
external rectenna areas of 0, 500, and 1000 ft 2 


15 


respectively. Since the design chosen has no external 
rectenna it must lie on the zero rectenna area line. 

Each series of graphs represents a different combination 
of an N lt (Ultimate load factor, the point where 
failure occurs) of 2, 2.5, or 3 with an aspect ratio of 

15 or 25. 

The benefit of raising the N ult is that the 
aircraft structure would be better able to withstand 
loads induced by gusts or maneuvering. At low 
altitudes, a small gust of wind can cause a significant 
increase in the wing loading resulting in deformation or 
failure of the wing. The tradeoff is that the wing 
structural weight and aircraft gross weight increase 
with the ultimate load factor. The effect of this 
increase is compounded by increased drag and rectenna 
area. An unexpected effect of reducing the aspect ratio 
is demonstrated by the graphs. The expected lessening 
of the lift-to-drag ratio and wing area did not 
materialize due to the unique relationship of the design 
parameters. Thus the reduction of the aspect ratio from 
25 to 15 does not produce a significant drag increase 
but results in a reduction of wing weight. 

4.3 Selection of Design 

A parametric study was conducted using the weight 
estimation code. All data was generated at a constant 
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wing thickness-to-chord ratio of 12%. The aspect ratio 
of the wing was varied from 25 to 15 while the aspect 
ratio of the horizontal tail was fixed at 25 and the 
vertical tail was fixed at 5. The ultimate load factor 
was varied between 2 and 5. Data obtained at higher 
ultimate load factors (greater than 3) converged only at 
unrealistically high design lift coefficients <C L > 1.0) 
which are not obtainable with a low Reynolds number 
airfoil. The optimum lift coefficient is .6 in level 

flight. 

The maximum acceptable wingspan is limited by the 
need to keep the wing tip deflections low. Therefore, 
the wingspan was constrained to a value less than 250 
feet. This limitation is based on the total wing weight 
allotment which is not enough for the heavy structural 
support needed to accomodate a larger wing span. It was 
also necessary to have the ultimate load factor as high 
as possible to maintain structural integrity and damage 
tolerance. Due to the increase of wingspan with aspect 
ratio, the final design reflects a compromise in design 
aspect ratio needed to accommodate the span limitations. 
The final configuration does not encompass any external 
rectenna because the parametric study failed to reveal 
any benefits for such a configuration. However, two 
drawbacks are evident. First, an external rectenna 
would produce a large drag penalty and secondly there 
are structural problems inherent in attaching it to the 
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aircraft. The wing area is all that is necessary to 
absorb the power required. The tail area can be used to 
absorb extra power that will be lost due to the angle of 
incidence between the beam and rectenna caused by the 
flight bank angle. The component weights for the first 
iteration are given in Figure 4.8. These values 
describe the configuration which has subsequently been 
analyzed in detail during this project. 

4.4 Comparison with Conventional Aircraft 

The major difference between this vehicle and a 
conventional aircraft is the direct relation between the 
wing area and available power. This connection results 
in parametric limits which are different than what had 
been anticipated. The most significant limiting factor 
on the design has turned out to be the wing span. Fuel 
constitutes a significant percentage of the weight of a 
conventional aircraft. There is a direct weight savings 
from eliminating liquid fuel which also results in a 
weight savings in the other aircraft components. 

The unique relationship between wing area and power 
available produces another important variation from a 
conventional vehicle with an air-breathing engine. The 
optimal Mach number for this aircraft is much lower than 
what was initially expected. Whereas an air-breathing 
aircraft flying at 100,000 ft would be expected to 
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operate most efficiently around Mach .7, our optimum 
Mach number falls around Mach .45. This is due to the 
fact that the wing area increases with cruising speed in 
order to provide sufficient rectenna area to absorb the 
required power. The wing area is much larger than for 
an air-breathing aircraft design of comparable weight. 

At such high altitudes, a large wing area is need to 
produce enough lift in the low density air at low flight 
speeds . 

4.5 Configuration 


The aircraft configuration was chosen after the 
preliminary weight estimate was obtained. Data was then 
available on expected cruising speed, altitude and 
flight duration. 

The cruising speed of the aircraft was determined 
as Mach .44 at an altitude of 100,000 feet. This allowed 
for a wing with no sweep which also produced a wing 
weight savings in comparison to the swept configuration 
required for operation at higher design Mach numbers. 

The low density of the air at design altitude leads to a 
large planform with low wing loading. In addition, 
typical gust velocities of up to 40 feet per second 
resulted in large normal accelerations (10 to 30 g) at 
low altitudes. 

The mission will require a payload consisting of 
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sensors that will sample ambient air. It is important 
that the air not be disturbed by the approach of the 
aircraft. 

4.5.1 General 

First, the possibility of an all-wing configuration 
was considered. An all-wing configuration represents 
the maximum lifting area possible for an aircraft of a 
given size. This was important because the rarified air 
at cruise altitude will reduce lift potential. Since 
the aspect ratio of the wing is 15 external bracing 
might be necessary. This issue is discussed in section 
6 . 0 . 

Figurss 4.9 — 4.11 illustrate design concepts with 
a canard and aft wing. Because of possible directional 
instability, the wing was relocated further forward and 
a twin fuselage configuration was used because this 
arrangement provides greater structural rigidity and 
internal volume. 

The tail configurations examined included butterfly 
and conventional horizontal— vertical tailplanes and 
canards. A canard was chosen because of its beneficial 
lift characteristics and an aft tail was included 
because of its effect on structural stiffness and pitch 
control at high altitude. The butterfly tail 
configuration was dropped when the use of a second 
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fuselage was made. 

The microwave propulsion system does not require 
fuel on board. The controls will be fly-by-wire, so the 
control surfaces will use electrical actuators, which 
means no hydraulics will be needed. The electrical 
system is representative of the wiring needed to take 
power from the rectennas to the engine. With a ground 
control system in use, the avionics will be grouped with 
the remote control system. The list of components 
includes the wing, horizontal and vertical tails, the 
fuselages, electrical wiring, payload, engines, 
propellers, controls and landing gear. 

The next components to be placed were the 
propellers. Three locations were considered. First, the 
placement of the propellers at the rear of the fuselages 
would be beneficial because the atmospheric sensors 
would be free of the turbulence created by the 
propellers. The drawback of this configuration is that 
propellers would be in the slipstream of the main wing 
and tail surfaces. This will result in propeller 
vibrations and a reduction in propeller efficiency. The 
second location would again be aft, but situated atop 
the vertical tail, free of the wing's wake. In this 
position the propellers are only affected by the airflow 
over the tail surfaces, providing a moderate 
improvement. This location is not considered feasible 
because of the offset thrust line which would cause a 
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large negative pitching moment. Another problem would 
be a complicated drive train because the power would 
have to be routed from the engine located in one of the 
fuselages through the vertical tail to the propeller at 
the top. The last choice was to locate the propellers in 
the front of the airplane. This allowed for a simple 
drive train. 

The landing gear was positioned fore and aft on 
each fuselage. Outriggers can be located on the wing 
tips if needed. A minimum ground clearance for the 
propellers of 2.25 feet was chosen, permitting a 
rotation angle of 15° as shown Figure 4.12. The final 
configuration is reproduced in Figure 4.13. 


4.5.2 Effects of Planform 

The planform of a wing is collectively comprised 

of the sweep of the leading edge, aspect ratio, AR, 

taper ratio and the general shape of the wing as shown 

in Figure 4.14. The selection of the planform 

significantly influences the vehicle aerodynamics. The 

aspect ratio affects the aircraft in two ways. An 

increased aspect ratio gives improved drag 

characteristics as demonstrated from the induced drag 

C » C 2 /rrAR. However, the wing weight is directly 
Di L 

proportional to the aspect ratio. This effect is 
illustrated in the weight estimation computer code which 
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shows a wing weight variation proportional to AR 
The taper ratio affects the vehicle aerodynamics by 
approximating an elliptical planform. Such a planform 
would yield a minimum induced drag. The wing sweep acts 
to increase the critical Mach number, which is the 
freestream Mach number where local velocities reach 
sonic speed. As the Mach number approaches M cr there is 
an associated rise in drag. 

4.5.4 Selection of Planform 

The design conditions for the microwave powered 
aircraft leave limited room for planform selection. For 
low speed flight there is little advantage to a low 
aspect ratio or high sweep angle. The operation of this 
vehicle fits into this category with a subsonic flight 
speed of Mach 0.44 at cruise altitude. An aspect ratio 
of 15 was selected based on the data in the weight 
estimation parametric study. This value represents a 
compromise between wing weight, aerodynamic efficiency 
and the limitation of the span length to 250 feet for 
structural reasons. The planform includes a slight 
taper on the outer 40% of the semi-span. Figure 4.15 
shows the actual planform selected. 

4.6 Summary of Aircraft Sizing and Configuration 
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After an initial analysis was performed by each 
design group based on the first iteration data, several 
questions arose. The weights from the first iteration 
seemed very low. The structures group felt the wing 
weight allocation was not enough to support a span of 
234 feet. The performance group desired more power for 
climb which would require a larger wing area. 

As the design evolved and the first iteration was 
completed, it became evident that a second iteration of 
the component weights was necessary. The new 
configuration contained a canard, two vertical tails and 
two booms instead of a single fuselage. A standard NACA 
airfoil of 8% thickness was chosen for the canard and 
horizontal tail and a NACA 10% thick airfoil was chosen 

for the two vertical tails. 

The weight estimation code was modified as follows 
to accommodate these changes: The aspect ratios and 

thicknesses of the canard, horizontal tail plane and 
vertical tails are based on the first iteration data. 

The values for these parameters are given in Figure 
4.16. An equation similar to the one for a horizontal 
tail was added to account for the weight of the canard. 
The fuselage equation weight was doubled to account for 
two booms each 50 ft. long and 4 ft. in diameter and the 
vertical tail equation weight was also doubled to 
account for the two tails. The design condition of C L 
.6 and the constraint of having no external rectenna 
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remain the same, however, the size of the main wing was 
not fixed due to the fact that it determines the power 

available to the aircraft. 

The second iteration deviates markedly from the 
initial estimate. The gross weight has more than 
doubled. This can probably be attributed to the 
increase in wing size caused by the increased weight and 
drag of the booms, canard and vertical tails. As 
mentioned previously, the wing weight equation is very 
sensitive at low values of wing loading. The overall 
wing span has also gone up considerably. In the first 
iteration, the span was fixed at 250 ft. to limit 
deflections. For this second iteration, the limit is 
relaxed due to structural findings that deflection are 
not as great as expected, and the possibility of 
reinforcing the wings with compression members. The 
second iteration component weights are illustrated in 
Figure 4.17. 
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5.0 Aerodynamics 


Using the results of the weight estimation program, 
it was found that the aircraft lift coefficient needed 
to range from 0.6 to 1.0, at a design Mach number of 
0.44. At low altitudes, the selection of airfoils which 
meet these requirements is simply accomplished by 
choosing a standard airfoil. At the design altitude of 
100,000 feet, however, due to the rarefied air, the 
chord Reynolds number is fairly low (about 711,000). 

At low Reynolds numbers, the transition of the 
boundary layer from laminar to turbulent typically 
involves a transition bubble. The transition bubble 
consists of a boundary layer separation just before 
transition, followed by reattachment and turbulent flow 
as shown in Figure 5.1.1. In the design of airfoils, it 
is generally desirable to keep the boundary layer 
laminar for as long as possible, as a laminar boundary 
layer induces significantly less skin friction drag than 
a turbulent one. However, a turbulent boundary layer 
has a much greater resistance to separation and may be 
advantageous at points where a large adverse pressure 
gradient exists. Therefore, it is desirable for the 
boundary layer to transition to turbulent at a point on 
the airfoil before a large adverse pressure gradient can 
induce laminar separation. The pressure distribution 
which the airfoil is designed to match, or the target 
pressure distribution, shown in Figure 5.1.2, features a 
flat upper surface in the transition region in order to 
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minimize laminar separation. 

In designing an airfoil for minimum drag, 
therefore, the target pressure distribution must be such 
that laminar flow is preserved as far aft as possible. 
Transition must then be induced before a steep adverse 
pressure gradient can cause flow separation. This can 
be accomplished by designing a pressure distribution 
with a flat rooftop and a long transition ramp before 
the major pressure recovery as shown in Figure 5.1.2. 

The transition ramp consists of a very gentle adverse 
pressure gradient which will induce transition at some 
point long its length without causing laminar boundary 
layer separation. In this way, the amount of laminar 
flow will be maximized, and skin friction drag 
minimized. 

The designed pressure distribution must be as 
smooth as possible. Otherwise, the resulting airfoil 
may have irregularities which cause premature transition 
of the boundary layer. This would result in an 
undesirable increase in the drag coefficient. 


5.1.1 Elliptical Wing Loading 


Ideally, an elliptical wing loading would be 
designed to provide for minimum induced drag. In order 
to acheive an elliptic wing loading with a non-elliptic 
planform, the airfoil section lift coefficient must be 
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The airfoil section lift 


varied along the span, 
coefficient is calulated based on the selected wing 
planform geometry from the classical elliptic lift 
equation given in Reference 13. 


o, =. yri - (2Y/'b ) ->i (5.D 

C(y) 

Where r is the circulation, y is the location along the 
semi-span which has a value of zero at the fuselage, 
is the free stream velocity, b is the overall span 
length and c is the local chord length which may be a 
function of location. 

This expression is evaluated by calculating r 
(circulation) from: 
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r = 


(5.2) 


/>» V . b » 

Where L is the lift force which equals weight in level 
flight and p „ is the air density. Figure 5.1.3 shows 
the airfoil section lift coefficient variation along the 
span. This lift distribution has an average value of .6 
over the entire span, which is the C L for this design. 
Using this figure, three airfoil sections were designed 
to be installed at specific locations along the span to 
meet the local airfoil section lift coefficient 
requirements . 
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5.1.2 Airfoil Design Overview 


The airfoil design procedure is primarily one of 
trial and error. There is no precise mathematical 
method from which a good target pressure distribution 
can be prescribed. The best tool in designing airfoils 
is an experienced designer and data collected by 
previous experiments and designs. To begin the airfoil 
design process, a baseline target pressure distribution 
was sketched. This target is based on data from 
References 10, 12, 19, and 23 on low Reynolds number 
airfoil design. 

Figure 5.1.4 shows the baseline airfoil section 
pressure distribution (L - Series Airfoils), which was 
then modified twice on the upper surface to produce a 
total of three variations. The first modification (P- 
Series Airfoils) features the rooftop further away from 
the leading edge, a more favorable pressure gradient on 
the upper surface at the leading edge and a shorter 
transition ramp. The second modification (Q-Series 
Airfoils) has the rooftop closer to the leading edge and 
a longer and steeper transition ramp. 

An airfoil section was generated for each of the 
three target pressure distributions using the computer 
code CPINVRS . FOR of reference 15. Figure 5.1.5 
describes the P-series designed airfoil. Figure 5.1.6 
gives a similar description for the L-series airfoil. 
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Finally, Figures 5.1.7 illustrates the Q-series designed 
airfoil. Each of these variations has a = 0.74, 
which is the section lift coefficient for the root of 
the wing. The L, P, and Q airfoils were then analyzed 
using the computer code Grumfoil of Reference 20. The 
results of this viscous compressible analysis show that 
increasing the transition ramp length will decrease 
drag. Also, a more gentle pressure gradient on the 
upper sufrace will delay laminar separation. Using this 
data, the L-series airfoils were chosen for maximum 
performance . 

Based on the spanwise lift distribution, three 
spanwise locations were selected and airfoil sections 
were designed for the appropriate lift coefficients by 
scaling the initial L-series target. The first section, 
LI is designed for the wing root. The next airfoil, L2 
is designed for installation at the juncture of the wing 
inner and outer panels. The final airfoil, L3 is 
designed for the wing tip. The design lift coefficients 
for these airfoils are 0.74, 0.59, and 0.40, 
respectivly . Figure 5.1.8 gives the LI, L2 and L3 
target pressure distributions. The LI, L2 and L3 
geometry is given in Figures 5.1.9-11. 

5.1.3 Design Procedure 

Each of the three L-series target pressure 
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distributions was first analyzed using thin airfoil 
theory to provide an initial estimate of the airfoil 
shape. The computer program DESIGN. FOR of Reference 13, 
based on thin airfoil theory, was used for this purpose. 
The input for DESIGN. FOR is in the form of the upper and 
lower surface target pressure coefficients specified at 
locations along the chord. The output airfoil 
coordinates are then used as a starting airfoil guess 
for the inverse program CPINVRS.FOR, Reference 15, a 
thick airfoil design code with uses the Theodorsen 
method to calculate pressures under the assumption of 
inviscid, incompressible flow. This program iteratively 
determines the required airfoil sections, providing 
airfoil coordinates at 25 chordwise locations on both 
the upper and lower surfaces. 

5.1.4 Airfoil Analysis 

The airfoils generated using CPINVRS were first 
analyzed using a computer code which provided an 
inviscid, incompressible analysis. Resulting pressure 
distributions are shown in Figures 5.1.9-11. In order to 
obtain drag data and more accurate pressure data, the 
airfoils were then analyzed using Grumfoil, a code which 
provided a viscid, compressible analysis of the 
coordinates. In Figure 5.1.12, the lift curves obtained 
from the inviscid analysis are shown. These are seen to 
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differ from those found in the viscous analysis, shown 
in Figure 5.1.13. The nonlinearities in Figure 5.1.13 
arise from viscosity and compressibility effects. 

Since the design procedure is based on an inviscid, 
compressible analysis, the airfoil sections were 
designed with zero thickness trailing edges. In some 
cases, as in Figure 5.1.14, this induces an abrupt 
change in surface slope adjacent to the trailing edge, 
which causes increased drag and unreliable convergence 
of a viscous, compressible analysis. To avoid this, the 
airfoil coordinates at points near the trailing edge 
were extrapolated as seen in Figure 5.1.14. This 
extrapolation created a trailing edge of finite 
thickness, which significantly reduced the drag in the 
viscous compressible analyses . 

When analyzed under free transition and at 
realistic Reynolds numbers, the L-series airfoils were 
found to have the best lift and drag characteristics at 
design conditions. This analysis was performed at the 
design Mach number of 0.44. Grumfoil analysis of 
sections LI, L2, and L3 shows that the incompressible 
analysis is fairly accurate, although there are some 
slight deviations in pressure at the leading edges. 
Design point pressure distributions for LI, L2, and L3 
are shown in Figures 5.1.15-17. Airfoil section drag, 
shown in Figure 5.1.18, is a significant improvement 
over the previously designed airfoils. 


33 


5.1.5 Summary of Airfoil Design and Analysis 

Comparison of the L-series airfoils against typical 
NACA airfoils of comparable thickness shows that the L 
designs have significantly lower drag near the design 
point. Drag polar comparisons are shown in Figures 
5.1.19-21. Also, the large amount of camber in the L- 
series airfoils results in a higher lift per degree 
angle of attack, as shown in Figures 5.1.22-24. These 
airfoils provide the necessary performance for the 
aircraft, although further refinement of the target 
pressure distribution could give even better airfoil 
section performance. 
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5.2 Wing-Body Analysis 


The objective of the wing-body analysis was to 
provide an aerodynamic study of the wing-body lift and 
drag characteristics. This study was accomplished by 
selecting a wing planform, integrating the designed 
airfoil sections into the wing and generating a 
geometric twist distribution. The wing pressure 
distributions are analyzed by using the program Flo - 
30m. The results approximate an elliptical wing loading 
and maintain a minimum induced drag. 

5.2.1 Twist Distribution 

Based on the three airfoil sections designed, a 
wing twist distribution can be determined. This 
distribution allows the loading to approach elliptical 
by varying the effective angle of attack along the span. 
Once this is done, each airfoil section will produce the 
lift coefficient described by Figure 5.2.1 For each of 
the three airfoils designed, the equation for the linear 
portion of the c . vs a curve is determined based on the 
results of the inviscid, incompressible analysis 
performed with CPDIST.FOR (Reference 14). The general 
equation for this portion of the graph is given in 
equation 5.3. 
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(5.3) 


c, = m a -- + c, 

1 eff l»-o 

Here, is the section lift coefficient, in is the lift 
curve slope, a eff is the effective angle of attack in 
radians and c. is the zero angle of attack section 

la-0 

lift coefficient. Figure 5.2.2 gives the lift curve 
equations for each of the L— Series airfoils. 

Since each airfoil is designed for a specific span 
location (LI at root, L2 at .6 semi-span, and L3 at .9 
semi-span) , the intermediate sections roust be 
approximated. As a first approximation, these sections 
gj-g evaluated as linear combinations of the two boundary 
airfoils. As the airfoils vary from root to tip the 
lift curve slope and the zero angle lift coefficient 
change. These two important parameters are also 
calculated using a linear approximation. The equations 
for the variation of these parameters are given in 
Figures 5.2.3 and 5.2.4. Substituting the relations 
from these figures into the general equation for c^ vs. 
a yielded a single algebraic equation that governs the 
lift curve over the span. Combining this equation with 
the section lift coefficient variation gave the 
effective angle of attack along the wing span as shown 
in Figure 5.2.5. The effective angle of attack near the 
root is not accurate since this analysis was performed 
without regard to fuselage interference effects. 

A first order linear approximation of the twist 
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angle holds the root section fixed at zero twist and 
varies the twist linearly in a nose down direction to 
the wing tip as shown in Figure 5.2.6. However, the 
effects of fuselage interference at the wing root 
necessitate increasing the root incidence angle to some 
positive (nose up) value. An elliptical distribution 
requires the tip section lift coefficient to be zero. 
Although the tip section is set at the zero lift angle 
of attack, according to the analysis, lift does remain 
on the tip. 

A second order linear approximation was determined 
to more closely follow the actual twist distribution as 
shown in Figure 5.2.7. Figure 5.2.8 shows the actual 
wing twist of the airfoil sections at several span 
locations. Figure 5.2.9 shows the airfoil twist 
variation viewed from the tip to the wing root. 

5.2.2 Compressible Wing-Body Analysis 

Using the wing geometry, airfoil sections and twist 
distribution described above, a compressible wing-body 
analysis was carried out using the program Flo-30m. The 
input to Flo-30m consists of airfoil sections and twist 
angles at various span locations. The computer code will 
interpolate to generate the remainder of the wing 
sections. After the wing is entered, the body geometry 
is described in similar fashion. Three iterations were 
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performed on the wing-body combination in level flight. 
The first was the untwisted wing geometry. The next two 
iterations were based on the first and second twist 
approximations. Only the data for the untwisted and 
final geometry is presented. 

The initial Flo-30m analysis was performed on an 
untwisted wing geometry. Flo-30m was run over a range 
of angles of attack from -2.0° to 6.0°. Using the C L 
versus a curve (Figure 5.2.10) the geometric angle of 
attack required to produce the design lift coefficient 
of .6 was determined (a=.375°). By analyzing the wing 
at this angle and examining the lift coefficient 
variation along the span a comparison can be made to the 
elliptical lift variation as illustrated in Figure 


5.2.11. 

Figure 5.2.11 illustrates the effects of fuselage 
interference close to the wing root. A large effective 
angle of attack is needed near the fuselage to approach 
the desired lift coefficient. The theoretical lift 
curve slope corrected for Mach number, sweep and aspect 
ratio for a finite wing was calculated from equation 5.4 
given on p 11—2 of Reference 24. 


dC, 


da 


m * AR 

2 + /t4 + AR 2 B* ( 1 + tan 2 0/B 2 ) ] 


(5.4) 


Where m is the airfoil lift curve slope corrected for 
finite span, AR the wing aspect ratio, B is a correction 
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factor for Mach number given as J { 1 - M J ) and 9 is the 
sweep of maximum thickness line. When the calculation 
was performed the corrected wing lift curve slope of 
6.99 per radian matched the Flo - 30m output. 

To improve the results, the wing was analyzed using 
the second approximation twist distribution. The 
versus a curve for the twisted geometry is also shown in 
Figure 5.2.10. The geometric angle of attack 
corresponding to the design lift coefficient was 
determined (a = . 2 °). The improvement in the lift 
distribution can be seen in Figure 5.2.11. The fuselage 
effects were still prominent. Figure 5.2.12 displays 
the wing loading for the elliptical, untwisted and 
twisted geometry. The twisted geometry more closely 
represents the elliptic wing loading and had a smaller 
induced drag. The lift curve slope for the twisted wing 
remained at 6.99 per radian. The twist however shifted 
the entire lift curve upward (Figure 5.2.10). 

The wing-body analysis was performed using level 
flight conditions. The actual flight path of the 
aircraft will be circular as described in section 2.0. 

In a circular flight path the effective wing lift 
coefficient required to maintain altitude will be 
slightly higher (see section 9.0). 

The Flo - 30m results for the untwisted wing in 
level flight are given in Figures 5.2.13 - 5.2.16. These 
figures show the pressure distribution on the wing at 
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various angles of attack. Similar diagrams for the 
twisted wing are shown in Figures 5.2.17 - 5.2.20. At 

the design lift coefficient, where the wing loading is 
close to elliptical, the pressure distributions are 
very favorable. At the off -design conditions, 
especially at high angles of attack, suction peaks 
develop on the wing tip. These suction peaks are due to 
the twist distribution being developed for specifically 
the level flight design C L and will cause the wing tip 

to stall. 

5.2.3 Theoretical Drag Calculations 


Following the methodology in Reference 24, 
theoretical calculations of the wing-body aerodynamics 
were carried out. Appendix A. 4 gives the details of 
these calculations. Figure 5.2.21 shows the theoretical 
wing-body drag polar. The C DMIN value from this pola 
was higher than the value used in the weight estimation 
program. This discrepancy lead to the revision of the 
drag polar in the weight estimation code for a second 
weight iteration. The viscous drag factor, K" , is much 
greater than the .007 value used in the initial 
analysis. A larger value of K" made the drag polar 
narrower. Therefore a small change in lift coefficient 
produced a large change in drag. This effect led to 
rapid drag increases. 
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Looking at the lift coefficient versus lift-to-drag 
ratio (Figure 5.2.22) demonstrated the actual L/D ratio 
is close to the weight estimation code L/D ratio. In 
fact, the plane would operate at maximum L/D at the 
level flight design lift coefficient of 0.6. Figure 
5.2.22 also confirms the rapid drag increases when 
deviating from design C L - The L/D ratio at C L = 0.6 is 
approximately 25. When in a circular flight path the 
L/D ratio decreases to 23.5. 


5.2.4 Summary of Wing-Body Analysis 


The wing-body analysis provided a twist 
distribution which approximates an elliptical wing 
loading. From the Flo - 30m analysis, it is clear that 
the low-Renolyds number wing design will produce enough 

lift for level flight (C^ = 0.6 at a g eo]in ® ^ ^ ’ 
on the theoretical drag calculations, the aircraft will 
have a L/D ratio high enough to ensure that the power 
available will be sufficient to overcome the drag force. 
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6.0 Structural Design & Analysis 



The primary objective of the structural analysis 
was to design the structural configuration and perform 
the static and dynamic analysis using the finite element 
software package ANSYS. The structural analyses 
consisted of three progressively complex models. The 
first two models embodied just the wing while the third 
model embodied the wing, tail, canard and boom. The 
wing aerodynamic loading distribution on model 1 varied 
spanwise and was constant chordwise while models 2 and 3 
varied spanwise and chordwise as shown in Figures 6.1.1 
& 6 . 1 . 2 . 

6.1 Static Analysis 


All three models were analyzed statically. The 
first model consisted of an aluminum skin with an 
internal structure of graphite - epoxy. It was 
comprised of 66 nodes, 60 beam elements and 60 plate 
elements. The chordwise pressure distribution for this 
run was constant as shown in Figure 6.1.1. For a more 
detailed discussion on how the first model was developed 

see Appendix A. 5. 

After running this model the following results were 
obtained. A 16.5 foot vertical deflection resulted, as 
shown in Figures 6.1.3 & 6.1.4. This deflection was 
small enough be analyzed using small deflection theory, 
as discussed in Appendix A. 5. 
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The average 


maximum shear stress was 2184 psi, as 
shown in Figures 6.1.5 - 6.1.8. This stress was below 
the ultimate shear strength of aluminum alloy 2014-T6, 
of 42 ksi . Note that most of the shear is carried in 
the bottom skin surface due to the upward pressure on 
the wing. 

The average maximum equivalent stress was 33,153 
psi, as shown in Figures 6.1.9 - 6.1.12. Again this 
stress is well below the ultimate tensile strength of 
aluminum 2014-T6 as well as graphite - epoxy which have 
strengths of 70 ksi and 69.9 ksi, respectively. Most of 
the bending stress is carried in the I beams as 
illustrated in Figures 6.1.10 - 6.1.11. However, small 
amounts of bending stress are carried by the skin 
surface, as seen in Figures 6.1.9 & 6.1.12. The skin 
surface where stress was maximum was located near the 
centerline . 

With this configuration the weight of the wings is 
2050 pounds. This is only 76 % of the allocated weight 
of 2704.6 pounds, allowing for further strengthening of 
the wing. The wing weight allowance was determined from 
the weight estimation analysis in Section 4.0. 

The first model was then altered to take into 
account the chordwise pressure distribution as discussed 
in Appendix A. 5. Running these changes with a gravity 
load produced better results, reducing the vertical 
deflection to 14 feet. However, the chordwise deflection 
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was 7.33 inches, which was too high for this model. 

This deflection was due to the low stiffening in the 

wing. 

Since the first model was too flexible, a second 
model was developed. The resolution of the second model 
was increased to 168 nodes, 180 beam elements and 180 
plate elements, to more accurately represent the 
structure. The entire structure was made out of 
graphite - epoxy. This material was chosen because it 
reduced wing weight while increaseing the strength of the 
wing. The geometry and wing aerodynamic loading 
distribution are presented in detail in Appendix A. 5 and 

Figures 6.1.13 — 6.1.16. 

Using the second model, four different analyses 

were run, including a gravity loading of lg, a pressure 
loading, a pressure and gravity loading, and an 
externally braced configuration with both gravity and 

pressure . 

The results of the gravity loading and the pressure 
loading analyses are presented in Appendix A. 5. 

For the combined pressure and gravity loading the 
average maximum equivalent stress was 36,955 psi, the 
maximum vertical displacement was 9.4 feet, and a 
maximum horizontal displacement was 5.3 inches, as shown 
in Figures 6.1.17-6.1.25. The horizontal deflections 
were reduced by 30% and the stresses were still quite 
reasonable while maintaining the wing weight allocated. 
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The last case for this model was for an externally 
braced configuration. The objective of this run was to 
see if a brace would significantly reduce the stresses 
and deflections. The brace was modeled using a one 
element bar with a modulus of elasticity of 20E10 psi to 
represent a rigid element. The external brace was 
connected to the middle I beam at a distance of 40.75 
feet out from the root and to the upper tail as shown in 
Figure 6.1.26. The average maximum eguivalent stress was 
lowered from 36,955 to 27,039 psi. The vertical 
deflection was reduced from 9.4 to 3 feet. The shear 
stress was increased to -8866 psi, approximately 5 times 
as much as without the brace. These results are shown in 
Figures 6.1.26-6.1.33. The compressive forces within 
the brace were 26.2 kips. These forces tended to 
produced buckling, as disscussed in Appendix A. 5. 

In order to construct a brace to withstand the 
given load, a larger cross - section was needed. The 
minimum acceptable diameter, using graphite - epoxy, was 
calculated to be 5 inches. However, this cross - section 
would generate unacceptable drag. This could be overcome 
by placing an airfoil around the braces. The airfoil 
would have a mean chord of about 50 inches. However, 
this would drastically alter the design of the aircraft. 

The third and final model was comprised of the 
second model plus a canard, tail, and boom. The 
pressures applied to the canard and tail are shown m 
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Figure 6.1.34. In this model a pressure with gravity 
loading was analyzed. The maximum vertical deflection 
was 8.5 feet, as shown in Figures 6.1.35 & 6.1.36. The 
maximum shear stress is -1588 psi as seen in Figures 
6.1.37 & 6.1.38. The maximum equivalent principal 
stress is 32,609 psi as shown in Figures 6.1.39 - 6.1.42. 

6.1.2 S umma ry of Static Structural Analysis 

All three models show that this aircraft at an 
altitude of 100,000 feet can structurally withstand 
flight conditions of 3 g's. However, at sea level, 
flight conditions were determined to cause a 24 g 
aerodynamic wing loading. No analysis was performed on 
a 24 g wing loading, however, one can readily see that 
the structure would fail before it ever reached such 
conditions. The only possible way the plane could 
operate at sea level conditions is if a gust alleviation 
system is incorporated in the aircraft, as described in 

Section 9.0. 
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6.2 Vibration Analysis 


Dynamic analysis was required to study the stresses 
under the unsteady loading conditions and to ensure that 
flutter is not encountered in the operating envelope of 
the vehicle. Dynamic analysis was used to determine the 
natural frequencies and mode shapes of the wing 
structure. It was important to review this analysis 
because the vibration , which can be caused by 
aerodynamic excitations, occurs during flight. The wing 
will deflect in vibrating motions and can cause a 
failure in the structural design. In the analysis, 
three configurations were analyzed by utilizing the 
computer code, ANSYS. The first wing structural model 
was the preliminary design, the second wing structural 
model was a revised design of the first model, and the 
third model included an external brace. 

Major parts of the dynamic analysis are the free 
vibration analysis of the wing as the clamped free beam 
and the free-free analysis. 

6.2.1 Clamped Free - Beam Analysis 

The half wing structural model was developed 
with ANSYS. A free vibration analysis was performed 
with the wing as a clamped free beam. The clamped free 
beam configuration has the same boundary condition as 
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a cantilever beam. Since the fixed end can not have any 
displacements, the wing root was fixed by restraining 
three translations and three rotations. The behavior of 
the ANSYS model was compared with the results of uniform 
cantilever beam. 

6.2.2 Preliminary Wing Structural Design 


The preliminary wing structural design and model 
consisted of two spars with ten spanwise rib sections, 
as shown in Figures 6.2.2 - 6.2.16 and as discussed in 
Section 6.1. The spars, which were modeled as the I- 

beams, were tapered along the span. 

Natural frequencies (Eigenvalue) of the each mode 

are shown in Figure 6.2.1 - 6.2.16. 

Figures 6.2.2 - 6.2.4 illustrate the first bending mode 
of the wing; Figures 6.2.5 - 6.2.7 display the first 
torsional mode; Figures 6.2.8 - 6.2.10 reveal the 
second bending mode; and Figures 6.2.11 - 6.2.15 show a 
combined mode involving bending and torsion. Note that 
there was no bending in the plane of the wing as shown 
in the Figures 6.2.4, 6.2.7, 6.2.10, 6.2.13, and 6.2.16. 

The first mode frequency of the preliminary structural 
design was 0.858065 Hz. This value was extremely low. 

It also has insufficient stiffness in the structure to 
prevent flutter. 
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6.2.3 Revised Wing Structural Design 


The revised wing structural design was developed 
from the preliminary model as discussed in Appendix A. 5. 

Natural frequencies of the each modes are shown in 
Figure 6.2.17 - 6.2.32. Figures 6.2.18 - 6.2.20 show 
the first bending mode of the wing; Figures 6.2.21 - 

6.2.23 illustrate the second bending mode; Figures 

6.2.24 - 6 . 2 . 2.26 display the first torsional mode; 
Figures 6.2.27 - 6.2.29 show the third bending mode; 
and Figures 6.2.30-6.2.32 display the second torsional 
mode. Again, note that there is no bending along 
longitudinal direction as shown in the Figures 6.2.20, 

6.2.23, 6.2.26, 6.2.29, and 6.2.32. 

Compared to the preliminary model, the additional 

structure gave more stiffness to the wing and resulted 
in higher frequency values for each of the respective 
modes. However, these frequencies were still quite 
indicating that the wing had excessive flexibility. To 
assess the adequacy of the wing structure with respect 
to flutter, it would be necessary to construct the power 
density spectrum of the in flight aerodynamic loadings 
and then conduct a dynamic analysis to calculate the 
magnitude of the structural response and associated 
stresses. This analysis is beyond the scope of the 
project, in which we sought merely to identify and 
investigate the wing f ree-vibration modes. 
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6.2.4 Revised Wing with External Brace 


The result of the revised wing structural design 
showed improved eigenvalues of modal frequencies. 

However, it was still not satisfactory. Therefore, an 
additional external brace was design to hold the wing 
and give additional stiffness to the wing structure. 

Such a brace was connected from the tailplane to the 
mid-section of the half-wing as shown in the Figures 

6.2.34 - 6.2.48. 

Figures 6.2.33 - 6.2.48 show the eigenvalues for 
each mode. Figures 6.2.34 - 6.2.36 show the first 
bending mode of the wing; Figures 6.2.37 - 6.2.39 
illustrated the second bending mode; Figures 6.2.40 
6.2.42 show the first torsional mode; Figures 6.2.43 

6.2.45 display the third bending mode; and Figures 

6.2.46 - 6.2.48 show the second torsional mode. 

Compared to the previous model, the natural 

frequency of the first mode was significantly increased 
from 0.9377 Hz to 1.5486 Hz. However, the brace did not 
adversely affect other modes because there were no 
significant increases in other frequencies. Figures 6.2.39 
and 6.2.45 clearly show that brace caused an increased 
in-plane deflection of the wing, which was not apparent 
in previous models. The external brace did not improved 
the stiffness of the wing structure, therefore, the 
brace was not feasible for additional support. 
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6.2.5 Summary of Vibration Analysis 


The three analyses indicated that the natural 
frequencies of the wing were extremely low and the values 
were unacceptable. This wing structural design will 
fail. The wing must be redesigned with more stiffness 
in the internal structure in order to increase the 
natural frequency. 
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7.0 Stability and Control 


In order for this aircraft to fly its mission of 
non-piloted flight to 100,000 feet, it roust be equipped 
with a carefully designed control system. The first 
step in designing a control system is to define and 
identify the flight characteristics. An aircraft of 
this nature does not need maneuverability, but it 
requires a high degree of reliability; its systems must 
be designed redundantly. The level of stability 
inherent in the vehicle is the basis of its flight 
characteristics; it must have adequate stability, yet 
it must also have a sufficient amount of control, 
especially for landing. 

7.1 Static Stability Anaysis 


Stability may be thought of as the tendency for 
the aircraft to return to equilibrium after being 
disturbed. Disturbances may be generated by the flight 
computer's control actions or by atmospheric phenomena 
such as wind gusts, turbulence, or thermal convection. 
Two flight conditions are necessary for the aircraft to 
fly its mission successfully: it must be able to 
achieve equilibrium flight, and it must have the 
capability to maneuver effectively over a wide range of 
flight speeds and up to the design altitude. 

In order for the aircraft to maintain steady 
uniform flight, the resultant force and moment about 
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the center of gravity (CG) must both be equal to zero. 
This state is known as the trim condition. If the 
forces and moments do not sum to zero, the airplane 
will be subjected to translational and rotational 
accelerations. To achieve static stability then, the 
vehicle must develop a restoring force and/or moment 
which tends to bring the vehicle back to equilibrium. 

Figure 7.1.1 is a plot of vs. a. When Airplane 
1 is subjected to a disturbance which pitches the nose 
up, a moment develops which continues to pitch the nose 
up further. Clearly, Airplane 1 does not possess 
static longitudinal stability, for which a negative 
value of pitching moment derivative (dC^/dot) is 
required. 

Both Airplanes 2 and 3 have negative pitching 
curve slopes, yet Airplane 2 cannot be trimmed at a 
positive angle of attack (a). Thus in addition to 
having static stability, the airplane must also possess 
trim capabilities of positive angles of attack. 

A computer program was written to carry out the 
static stability computations. Each component's 
contribution to the pitching moment curves may be seen 
in Figure 7.1.2. It is noted that the fuselage and 
canard are destabilizing, while the wing and tail are 
stabilizing. The airplane pitching moment is stable; 
however, it will not trim at zero angle of attack in 
the absence of canard or tail deflection. 
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Using the data of Figure 7.1.3 and the analysis 
program (see Appendix A.7), the neutral point of the 
airplane was found to be 4.59 ft. aft of the leading 
edge of the wing. It was then determined that 
placement of the CG at 0.1 chords ahead of the neutral 
point would provide an acceptable level of static 
stability while allowing for adequate control. Thus 
the CG would lie 3.05 ft. aft of the wing's leading 
edge. Because there is no fuel carried aboard this 
aircraft, the CG position should remain constant during 
all phases of flight. Thus, at least in a preliminary 
design, it is not necessary to analyze the effect of 
center of gravity position on control. 

7.1.1 S umm ary of Static Stability 

The final task undertaken in the area of static 
stability analysis was that of generating pitching 
moment curves at different tail and canard incidences. 
These curves (Figures 7. 1.4-7. 1.9) define the limits of 
control of the airplane. 

At high angles of attack, an upward deflection of 
the canard will tend to stall the canard, creating a 
higher overall stability and a pitching moment towards 
equilibrium. In the regions where the canard stalls 
before the tail, the total canard contribution to 
stability disappears. This is seen as a decrease in 
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the slope of the airplane pitching moment curve. If the 
tail stalls first, its contribution likewise 
disappears — decreasing overall stability. This is seen 
as an increase in the slope of the pitching moment 
curve . 

The model for this analysis assumes that the lift 
curve slopes for the canard and tail are constant up to 
8°, beyond which the section lift curve slopes become 
zero and the lift coefficient is constant. (See data 
for NACA 64208 Airfoil, Reference 01). A more rigorous 
analysis requires a wind tunnel model to test for these 
nonlinear effects. 

Three types of horizontal control surface 
scheduling were examined; these correspond to Figures 
7. 1.4-7. 1.9. The first shows deflections of the tail 
while the canard is fixed at zero incidence. The 
second shows deflections of the canard while the tail 
is fixed at zero incidence. The third shows scheduling 
of equal deflection between the tail and canard. 

Another scenario for longitudinal control scheduling 
might be to deflect both surfaces in such a way that 
each donates an equal amount to the pitching moment. 
More elaborate methods might include surface 
deflections which are a function of the stresses in the 
structure, or of the flight dynamics. 

For the three scenarios examined, it was 
determined that the aircraft can trim over a range of 
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angle of attack from (-7.8° to 5.7'). This range 
corresponded to lift coefficients varying from (-0.58 
to 1.02). Figures 7 . 1.4-7. 1.9 demonstrate that this 
vehicle has adequate control regardless of the method 
of surface deflection scheduling, although it appears 
desirable to increase the maximum positive angle of 
attack at which the vehicle can be trimmed. 
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7.2 Dynamic Stability 


The purpose of this analysis is to determine 
whether or not the aircraft will be stable at all of its 
flight conditions. This was done by finding its natural 
modes of oscillation at given altitudes and checking 
whether they tend to converge (stable) or diverge 
(unstable) . 

The dynamic stability analysis for our aircraft was 
performed for altitudes of 100,000 feet and 25,000 feet. 
Cruise velocities are Mach .44 and Mach .09 at the 
respective altitudes . 

Only a longitudinal dynamic analysis was done for 
this aircraft, so the relevant factors were u (the 
dimentionless velocity change along the x axis), w (the 
dimentionless velocity change along the z axis), and q 
(the rate of change in pitch angle 0). 

The manifestation of an aircrafts response to 
disturbance occurs in two modes. The short period is the 
immediate response, and the long period, or phugoid, is 

the "long term" response. 

The first step in the analysis was the 
identification of the longitudinal stability 
derivatives. The derivatives were derived from basic 
equations of motion and small disturbance theory for a 
wing and tail (Reference 16). They were then modified 
to account for the aircraft's canard. 
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Certain assumptions were made concerning the 
derivatives. Compressibility effects were neglected. 
Direct contributions of thrust were neglected. 

The stability derivatives were calculated as shown 
in Figure 7.2.1. They were then put into the homogenous 
dynamic stability equation shown below in matrix form. 

(7.2.1) 


-du- 

dw 

dq 

kie J 


*u 

L «u +M w M z 




$>w z w 



-du- 

dw 

dq 

•-de- 1 
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«6 +M w Z 6 



/ * 
/ *T 


*The second matrix term is 0 since the analysis did 


not include control inputs 


Before the exact analysis was carried out, 
theoretical approximations for the long and short period 
frequencies and damping ratios were calculated using 
approximation equations derived from equation 7.2.1 
These theoretical values were then compared with the 
results of the exact analysis. 

Three types of analyses were carried out. The short 
period was calculated using a two-degree-of-freedom 
system of equations while the long period motion was 
found by writing a second order differential equation 
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for the velocity perturbation u. These equations are 
shown below. 


Short period equations: 
du = 0 

dw = (Z )dw + udq 

dq = (M+MZ )dw + (M +M w u 0 )dq 


} 


(7.2.2) 


Long period equation: 

u = (X u )u + (gZ u /u Q )u (7.2.3) 

The exact analysis was performed by solving 
equation 7.2.1 using the differential equation solving 
computer program ASDEQ. The main block diagram used as 
input for ASDEQ is shown in Figure 7.2.2. 

In order to find the natural modes of the aircraft, 
the control matrix was assumed to equal zero. This 
allowed the response of the aircraft to be calculated 
simply by changing the initial values of u, w, and q to 
disturb it from equilibrium. This simulates an induced 
perturbation in forward, vertical, and pitching motion 
respectively. The equations are then integrated with 
respect to time and the aircraft response over a chosen 
time period is shown graphically. The results for the 
airplane are given in Figures 7.2.3 - 7.2.5 and Figures 
7.2.6 - 7.2.8 for the examined altitudes ( 100,000 feet 
and 25,000 feet respectively). 
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7.2.1 S umm ary of Dynamic Stability 


The initial condition for the long period mode at 
both altitudes was u=2. The initial perturbation for 
the short period mode at both altitudes was a=l degree 
(.017 radians), as was that for the combined mode 
analyses . 

The separate long and short period results were 
generally predictable. The damping at 25,000 ft was 
much greater than that at 100,000 ft due to the greater 
dynamic pressure at the lower altitude. The airplane 
was stable at both altitudes, but the long period mode 
at 100,000 ft is just barely so. The short period mode 
at 25,000 ft can be seen to be more than critically 
damped in Figure 7.2.6 (it does not have an oscillatory 
response), while the short period at 100,000 feet can be 
seen in Figure 7.2.3 to be about 4 seconds. The 
difference in the long periods at the two altitudes is 
shown in Figures 7.2.4 & 7.2.7 to be 60 seconds at 
100,000 feet versus 12.5 seconds at 25,000 feet. 

The results of the combined mode analysis are far 
more surprising. The airplane is stable at 100,000 
feet, but slightly unstable at 25,000 feet (as shown in 
Figures 7.2.5 and 7.2.8 respectivly) . This is the 
reverse of what one would expect based on the relative 
dynamic pressures of the two flight conditions, as well 
as the opposite of what was predicted by the results of 


62 


the separate long and short period inodes. Regardless, 
the program was functioning correctly (Appendix A. 8) and 
the response derivatives were carefully checked, as was 
their precise input into the program. While an error is 
always a possibility, we consider it unlikely that this 
was the cause of the unexpected results. 

A possible explanation of the discrepancy is the 
extreme sensitivity of the aircraft to changes in 
apparent velocity at low altitudes. This may overcome 
the aircraft's tendency to recover from a disturbance. 
The sensitivity may be caused by an extremely low wing 
loading and may occur only through coupling of the short 

and long period modes. 
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8.0 Weight and Balance 



8.1 Component Sizes and Locations 


The weight estimation program provided weights of 
the components based on statistical data. Certain 
components such as the wing and fuselages had their 
weights revised during the structural analysis. The 
weight of the rectenna was a function of its surface 
area. This weight was estimated with reasonable accuracy 
in the weight estimation program. Other components such 
as the engines and gearboxes needed a second weight 
estimate based on a more accurate data. 

Each propeller was calculated to be 10^ feet in 
diameter using momentum theory and an assumed efficiency 
of 85%. Using this diameter and other data, the gearbox 
weight was calculated at approximately 100 pounds (see 
Appendix A. 6). The electrical motor was not designed, 
rather its physical dimensions were based on an 

existing commercial design. 

After the components were weighted, they were 
positioned in the aircraft in order to position the 
aircraft's center of gravity at a predetermined 
location . 

8.2 Mass Properties 

With the components sized and positioned, the 
mass moments and products of inertia about the 
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aircraft's center of gravity were calculated by 
assigning geometric shapes (flat plates, cylinders and 
flat discs) to the components and then using 
corresponding equations. Figure 8.1 is a side view of 
the aircraft with the components located and dimensions 
shown. Figure 8.2 is a data table showing all the mass 
properties of each component in addition to the X and Y 
coordinates (identified as X-Bar and Y-Bar respectively) 
of the aircraft's center of gravity. 


8.3 Summary 

The vehicle components were sized and located about 
the aircraft so that the center of gravity was 36.71 
feet behind the propellers along the roll axis of the 
aircraft. I was calculated to be 5.43E+5 slug-feet 2 ; 

X X 

I to be 4.41E+4 slug-feet 2 and I to be 517E+5 

yy 

slug. feet 2 by simplifying the component shapes into flat 
plates, cylinders and flat discs. 
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9.0 Performance 


9 . 1 Objectives 


The performance of an aircraft involves the method 
of takeoff, climb, cruise, and the path followed. The 
mission requirements of this aircraft require 
adjustment to substantial density variations resulting 
in a aerodynamic and structural tradeoffs. This 
section will analyze and discuss the intended flight 
conditions and parameters for the mission. 


9.1.1 Power 

The ground location of the power source will be 
fixed. Therefore, the aircraft must fly in a circular 
path at altitude about the power source. The power 
required and the availability of power at design 
altitude and also during the aircraft's climb to its 
cruising altitude must be analyzed in order to 
determine whether sufficient power is available to 

execute the desired flight path. 

The aircraft will not be flying in level flight 
but at a small bank angle, <p, in order to turn and 
remain near the ground power source. The power is 
transmitted most efficiently if the microwave beam is 
perpendicular to the rectenna. Flight in a turn with a 
stationary power source does not allow for this. The 
angle between the normal to the rectenna and the beam 
shall be referred to as the angle of incidence, 6, see 
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Figure 9.2.8. The availability of power will be 
decreased due to this angle of incidence in turning 

flight. 

The total power which must be transmitted to the 
aircraft is another performance consideration. The 
power source beams power in a conically shaped beam in 
which the density of the power decreases with distance 
from the center of the beam, as shown in Figure 9.2.19. 
The amount of power and the diameter of the transmitter 
must also be determined for the design conditions. 

9.1.2 Flight Path 

The flight path to the cruising altitude will be a 
spiral with the bank angle, flight radius, and velocity 
varying during the climb. A description of the path 
can be determined from the chosen flight conditions at 

each altitude. 

9.1.3 Climb 

The rate of climb, dh/dt, is a function of the 
available and required power and the aircraft weight, 

W, such that : 

dh/dt = ( Pavailable - Prequired ) / W (9.1) 

The amount of excess power available determines an 
aircraft's capability to climb. Since this aircraft 
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receives microwave power most efficiently in level 
flight, i.e. when the beam is perpendicular to the 
rectenna, an analysis of the excess power and climb 
performance was first constructed for a wings-level 
flight condition. 

A second analysis of the climb rate will consider 
flight in a turn. The corresponding bank angle, <p, 
angle of incidence, B, Mach Number, M, and flight 
radius, R, will be determined at each altitude after 
the maximum rate of climb the aircraft is capable of 
achieving is chosen. This will allow the optimum 
flight path to be described. 

9.1.4 Effect of Wind Gust 

The magnitude of the load factor, n, encountered 
due to a gust was determined in order to design the 
aircraft structure. A low wing loading, W/S, increases 
the load factor for a given gust velocity, w : 

n =l+ (p*a*w*V) / 2 (W/S) (9.2) 

Since the design condition of the aircraft at a 
altitude of 100,000 feet required a low wing loading, 
this can not be avoided. The lower air density at 
higher altitudes decreases the load factor experienced 
at a given airspeed. Encountering a gust at the design 
altitude will not be a performance problem since the 
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structure has been designed for a load factor of 3. 
However at lower altitudes, expected gust velocities 
produce very large load factors, far in excess of the 
values to which the structure has been designed. Thus 
a method of alleviating the effects wind gusts during 
the climb must be a part of the design. 


9 . 2 Power Analysis 


9.2.1 Level Flight 

In order to analyze the availability of power at 
the design altitude and during climb a computer program 
was written to calculate various characteristics of the 
aircraft as the speed and corresponding Mach Number 
increase. A listing of the code can be found in the 
Appendix A. 9. The outcome of the preliminary weight 
estimation were input into the program and the 
resulting coefficient of lift, lift, coefficient of 
drag, drag, thrust, propeller efficiency, power 
required, and power available were calculated as the 

Mach number and altitude varied. 

The power available and required at sea level, 
25,000 feet, 50,000 feet, 75,000 feet and 100,000 feet 
are shown in Figures 9.2.1 - 9.2.7. The results of 
this analysis indicate that there is sufficient power 
available at each altitude to maintain level flight. 
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9.2.2 Turning Flight 

For this aircraft, flight in a turn results in a 
reduction of available power and an increase in 
required power. The wings must produce more lift 

since : 


L = W / COS«p (9.3) 

In level flight the lift is equivalent to the weight, 
as the bank angle increases the necessary lift 
increases. The propeller efficiency, Ep, a function of 
the thrust, velocity, air density, and cross-sectional 
propeller area, also is reduced in turning flight: 

E p = 2/(1 + (1 + T/ ( 0 . 5 * p * V2 * Apropeller ) )l/2) (9.4) 

Power Available is reduced in a turn because the 
microwave beam is no longer perpendicular to the 
rectenna surface. The angle of incidence, B, is a 
function of the bank angle, <p, altitude, h, and flight 
radius, R (see Figures 9 . 2.8 -9.2.1) : 

B = v + 0 (tan© = R/h) (9.5) 

The power received by the motor is reduced by the angle 
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of incidence as follows : 


Preceived = Pd * S * cosB (9.6) 

The bank angle and angle of incidence which result in 
the maximum available excess power will be the optimal 
flight conditions, as long as the corresponding 
coefficient of lift is attainable. Since the speed and 
air density vary with altitude, so do the bank angle, 
angle of incidence, and radius of turn. 

The level flight performance program was altered 
to account for the changes taking place in turning 
flight. A power analysis of turning flight at bank 
angles ranging from 5° to 30° at five altitudes ranging 
from sea level to 100,000 feet was done using the 
revised program. Figures 9.2.11 - 9.2.18 show 
the power available and required in turning flight as 

the bank angle varies. 

The results of this analysis show that the 
aircraft has sufficient power to fly from sea level to 
100,000 feet in turning flight, although there is very 
little excess power at the design altitude. Additional 
power could be received if the canard and tail are used 
as rectenna surfaces. The optimal flight speed and 
radius increase with altitude as shown in Equation 

9.7 : 
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r = V2 / (g * tarup) (9.7) 


The coefficient of lift necessary to maintain flight 
also increases with altitude since it is inversely 
proportional to the air density . 

CL = L / (1/2 * p * S * V2 ) (9.8) 

At the lower altitudes, 5 to 25 thousand feet, the 
optimal flight speed is very low, in the range of Mach 
0.1. The necessary coefficient of lift is in the range 
of 0.4 to 0.5. As the flight altitude increases the 
necessary coefficient of lift also increases as well as 
the flight speed. It is difficult to design one 
aircraft that is capable of flight at sea level to 
100,000 feet since the air density varies considerably. 
A configuration having variable wing area is a possible 
solution to this problem. 

9.2.3 Total Transmitted Power 

The ground power source must be sized to 
accomodate the power requirements. The area of the 
transmitter and the total power transmitted are 
inversely proportional and a function of the power 
density (PDENSITY ), wavelength of the microwave 
radiation, altitude (h), and an efficiency factor (n) 

(Reference 6) : 
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At Pt = ( PDENSITY * wavelength2 * h2 ) / n (9.9) 


The density of the power has a parabolic power 
distribution in which the center has maximum power, see 
Figure 9.2.20. At the h beam width angle, BW1/2, there 
is only half of the power density. The maximum 
attainable power density is 700 W/m2. The calculation 
of total power for a range of transmitter areas is 
presented in Figure 9.2.21. 


9 . 3 Climb Performance 


9.3.1 Climb Rate 

The performance program was revised once again so 
that the climb rate could be calculated at the various 
altitudes. First the wings level climb was calculated 
and then turning flight, as shown in Figure 9.3.1. 

The amount of excess power versus the Mach number, 
shown in Figures 9.3.2 - 9.3.6, was used to determine 
the optimal flight speeds for each altitude. It was 
determined that the optimal bank angle was in the range 
of 5 C to 15° and that it varied with altitude. 

The program was run once again for bank angles 
ranging from 3° to 13° and altitudes ranging from 5,000 
feet to 100,000 feet. This information determined the 
climb rate and corresponding bank angle and flight 
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speed with the restriction of a CL < 0.8 at each 
altitude. Climb rate, bank angle, velocity, Mach 
number, coefficient of lift, and flight radius versus 
altitude are shown in Figures 9.4.2 - 9.4.7. 

9.3.2 Time to Climb 

The time to climb can be calculated using the 
inverse of the maximum rate of climb, dh/dt. This is 
calculated as : 


t = ffi dh (9.10) 

Jo (dh/dt) 

The minimum time to climb for both wings level 
flight and turning flight have been calculated and are 
given in Figures 9.3.7 and 9.3.8. 

9.4 Flight Path 


The aircraft must be towed to an altitude that it 
is capable of receiving the power from the ground power 
source without much of a reduction due to the angle of 
incidence. Twenty thousand feet is a sufficient 
altitude because there is an incidence angle of 
approximately 13° which does not result in a sizable 
decrease in power received by the rectenna . 

The flight path to the design altitude will be 


76 


determined from the speeds and bank angles resulting in 
the optimal rate of climb for each altitude, see 
Figures 9.4.1 - 9.4.3 and 9.4.5. The aircraft's bank 
angle will first decrease from 7 degrees at an altitude 
of about 20,000 feet to 5 degrees at an altitude of 

25.000 feet, stay fairly constant until an altitude of 

50.000 feet, and then increase steadily to 9 degrees at 

100.000 feet. The speed will be increasing during the 
climb from a Mach number slightly below 0.1 to a 
cruising Mach number close to 0.4 at 100,000 feet. The 
increasing speed accounts for the increasing radius, 
which varies from approximately 2,500 feet to 28,000 
feet at the design altitude. The result is a spiral 
flight path to an altitude of 100,000 feet and cruising 
flight in a turn at a bank of 9 degrees and flight 
radius greater than 28,500 feet. See Figures 9.4.1 
g # 4 # 7 to visualize the mission flight path. 

The number of turns the aircraft must make to 
climb to the design altitude can now be determined. The 
velocity at a given altitude is a function of the 
radial velocity and the flight radius : 

V = O * R (9.11) 

And the time to fly one turn, T, is a function of the 
radial velocity: 
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T = (2 * »)/n = (2 * W * R)/V (9.12) 


The change in altitude, d(alt), for one turn is 
proportional to the climb rate, dh/dt, at the altitude 

of flight : 


dh/dt = d(alt)/T (9.13) 
d(alt) = dh/dt * T (9.14) 

The aircraft will make 2 turns from 20,000 to 

25.000 feet. It will then continue to 50,000 feet 
making 8 more turns. The climb to 75,000 feet will 
require almost 8 more turns. The climb rate from 

75.000 to 100,000 feet is quite low, this increases the 
number of turns to 28. See Figure 9.4.1 to view the 
change in flight radius with altitude. The number of 
turns necessary is given in Figure 9.4.8. 


9 . 5 V-n Diagram 

A program was developed to calculate the V n 
Diagram which shows the load factor of the aircraft 
configuration due to a vertical wind gust. At five 
altitudes ranging from sea level to 100,000 feet the 
load factor was calculated for gusts ranging from -40 
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ft/s to 40 ft/s as shown in Figures 9.5.1 through 
9.5.5. The CL will never be greater than 1.5 or less 
than -1, therefore the associated load factors versus 
Mach Nu mb er for these CL's were calculated and plotted, 

as shown in Figures 9.5.6 - 9.5.10. 

The aircraft was structurally designed for a load 
factor of 3. If it had been designed for a load much 
greater than 3 the structural weight would have 
increased significantly. At altitudes above 60,000 
feet the present aircraft design is capable of 
encountering gusts above 40 ft/s. At lower altitudes, 
gusts as small as 10 ft/s result in load factors well 
above 3, see Figures 9.5.11 - 9.5.13. The flight 
speed at these altitudes is very low, yet the load 
factors induced by typical gusts will induce loads that 
exceed the design limitations of the wing. 

The aircraft is structurally capable of 
encountering very large gusts at high altitudes. This 
analysis indicates that as presently designed, the 
aircraft would require a gust load alleviation system 
at altitudes below 60,000 feet. This would take the 
form of sensors mounted on the aircraft which would 
activate the flight controls in response to a gust in 
order to prevent the development of load factors in 
excess of the design value. Alternatively, a variable 
wing area configuration might be a possible solution. 
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9 . 6 Results 


The results of this analysis showed that there is 
sufficient power available to climb in a spiral flight 
path and maintain turning flight at the design 
altitude . 

The coefficient of lift necessary to maintain 
turning flight at the design altitude is approximately 
0.8. The flight speeds at lower altitudes are quite 
low. At altitudes ranging from sea level to 50,000 
feet, a coefficient of lift ranging from 0.4 to 0.5 is 

adequate . 

The aircraft shall be towed to an altitude of 
20,000 feet by another aircraft. At 20,000 feet the 
aircraft will begin turning flight at a flight radius 
of approximately 2,500 feet using microwave power on 
the design spiral flight path. A very sensitive gust 
load alleviation system will be in use as the aircraft 
is towed and during the climb to 60,000 feet. Once the 
aircraft reaches 60,000 feet the gust load alleviation 
system will still be available but the sensitivity can 
be decreased. The aircraft will climb to the cruising 
altitude steadily increasing the flight radius until a 
altitude of 100,000 feet and radius of approximately 
28,500 feet is reached. This climb to 100,000 feet is 

estimated to take 3*75 hours* 

Approximately 15 to 25 MW of microwave power must 
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be transmitted to the aircraft from the ground power 
source. This amount of power is substantial and 

costly . 


9.7 Summary of Performance 

There are critical issuses in order to accomplish 
such a mission which our study has not yet resolved. 
Using the wing as a rectenna surface, the optimum wing 
loading is lower than would be required with some form 
of liquid fuel propulsion system. This low wing 
loading makes it difficult to design a structure which 
can withstand the large load factors induced by gusts 
at low altitude. Therefore, either a gust alleviation 
system or a variable wing area design should be 
considered . 

The total power required from the transmitter is a 
substantial fraction of a power plant. This amount of 
power is quite expensive and damaging to the 
environment. The cost and environmental impact will be 
discussed in Sections 12.0 and 11.0 respectively. 
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10.0 Beam Tracking System 



10.1 Overview 

There are two major engineering tasks associated with a 
microwave powered high altitude aircraft: the actual 
aircraft itself and the system that delivers microwave 
energy to the aircraft. The microwave energy could be 
transmitted to the aircraft from several locations including 
the ground, an aircraft at a lower altitude or even perhaps 
a satellite. The most feasible approach from an engineering 
standpoint is transmission from the ground. 

The transmitted beam of microwave power must be very 
narrow and accurate to maximize efficiency and minimize 
interference with other radio frequency sensitive objects. 
Since the microwave beam is very narrow, it must be kept 
pointing at the aircraft at all times. This is a 
significant engineering task when the aircraft operates at 

an altitude in excess of 70,000 feet. 

This chapter addresses the task of tracking a simulated 
aircraft with a beam of energy, by using a beam of ordinary 
light to simulate the microwave beam, a tracking system 
model undergoes design, implementation and evaluation from a 
controls engineering perspective and for its applicability 
to the "real world" problem. 

10 . 2 Background 

The concept of transmitting electrical power in a 
"wireless" fashion using microwave energy originated with 
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the early work on radio waves in the late nineteenth 
century. Heinrich Hertz and then later Nikola Tesla showed 
that energy could be transmitted through the air using radio 
waves. The work done by Hertz and Tesla utilized low 
frequency radio waves and proved that efficient power 
transmission would require shorter wavelengths. Radio waves 
with wavelengths less than a meter (frequencies above 300 
megahertz), known as microwaves, are required for proper 
beam focusing. The technology of the early 20th century 
did not allow generation of sufficient power levels in the 
microwave range for workable transmission of power. Major 
breakthroughs occurred during World War II when the klystron 
tube, and subsequently, the cavity magnetron were invented 
for use in radar systems. Both of these devices provided 
the necessary levels of power for feasible energy 
transmission using microwaves. 

In 1959, a team led by William C. Brown at Raytheon 
Corporation undertook one of the first practical 
applications of microwave power transmission. A small 
helicopter with a brushless DC motor became the first 
airborne vehicle powered by energy from a transmitted 
microwave beam. The helicopter utilized an array of small 
dipole antennas, each one connected to a diode rectifying 
circuit, which collected the microwave energy and converted 
it to DC. This type of rectifying antenna is known as a 
rectenna. The original helicopter, demonstrated in 1964, 
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flew to a height of 60 feet guided by tethers. Raytheon 
also fabricated a free-flying model which they successfully 
demonstrated to the Air Force, however the existing 
technology did not allow any practical applications. The 
most significant result of the program was the increased 
knowledge of energy collection using the rectenna, as 
discussed in Reference 5. 

The next major proposed application using microwave 
power transmission arose in the early seventies as a result 
of the energy crisis. The concept consisted of a network 
of satellites collecting solar energy, converting the energy 
into microwaves, then beaming the energy to the ground. The 
Solar-Powered Satellite (SPS) program, as it was called, 
proposed the use of satellites with collector panels the 
size of Manhattan Island. They would have generated up to 
300 billion watts of energy at a cost well into the 
trillions of dollars. Although the cost of the SPS program 
was prohibitive, approximately ten years of research 
provided significant advancement in the area of microwave 
power transmission. An improvement in rectenna technology 
emerged as one of the most significant achievements of the 
SPS research (Reference 5) . 

In an effort to maximize the efficiency of a microwave 
transmission system, William Brown (at Raytheon) and James 
Triner of NASA developed a thin-film rectenna from printed 
circuit technology. The thin-film rectenna weighs one-tenth 
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of the original rectenna developed at Raytheon and provides 
efficiencies in excess of 80 percent. The thin-film 
rectenna opened the door for aerospace applications of 
microwave power transmission, including the most recent 
breakthrough, the first microwave powered airplane 
(Reference 6) . 


10 . 3 Microwave Powered High Altitude Aircraft 

The first official flight made by an airplane powered 
by beamed microwave energy occurred on October 6, 1987. 

The airplane, called the SHARP-5 (Stationary High Altitude 
Relay Platform) , was designed and built by the University 
of Toronto Institute of Aerospace Studies (UTIAS) . Major 
funding from the project came from the Canadian 
Communications Research Center (CRC) because of their 
interest in the SHARP for communications applications 
(Reference 29) . 

The SHARP-5 aircraft has a fifteen foot wing span and 
carries a 3.5 foot circular disk aft of the wings for 
microwave reception. The entire underside of the aircraft 
is covered with thin-film rectenna made with diodes 
fabricated from silicon. The original thin-film rectenna 
utilizes diodes made of gallium arsenide which are more 
efficient for power conversion, but cost substantially 
more. 
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Power for the aircraft is generated by a battery until 
it reaches an altitude where it can intercept the microwave 
beam. The transmitted frequency of the microwave beam used 
is 2.45 gigahertz. There are two rectenna layers at right 
angles to each other to allow maximum energy collection 
while the airplane follows a circular path. 

An aerodynamical ly improved model. The SHARP-6B, has 
also been developed by UTIAS. The SHARP-6B is intended to 
be the prototype for a future full-scale version with a 130 
foot wing span and a 30 foot rectenna disk with 10,000 
dipoles capable of flying at altitudes above 50,000 feet at 

speeds near 120 mph* 

Future plans for autonomous aircraft include many 
applications. Radio and television broadcasts could be 
transmitted to the aircraft and then relayed down to a 
specific area, which would reduce the size of receiving 
dishes to a few feet. Possible uses in communications 
include telephone networking and cellular telephone service. 
Military reconnaissance and early warning of low flying 
aircraft and cruise missiles comprise other applications. 
Various type of airborne surveillance such as coastline 
monitoring could also be performed. The lower costs 
associated with a remotely piloted aircraft might make it an 
attractive alternative to satellites in certain applications 
(Reference 29) . 
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10.4 Microwave Transmission Systam 

The feasibility of a high altitude unmanned airplane 
powered by microwave energy depends largely on the cost of 
developing an effective transmission system. Current 
technology places the cost of a ground transmission system 
in the neighborhood of 30 to 50 million dollars, which is 
equivalent to the cost of launching a communications 
satellite. A microwave transmission system consists of the 
actual power transmitting elements and some control system 
that keeps the transmitting elements pointed at the target. 

The most popular transmission system concept currently 
pursued by both NASA and CRC involves the use of several 
thousand individual transmitting modules powered by standard 
microwave oven magnetrons. The appeal of this concept 
exists because of the availability and low cost of oven 
magnetrons, and the reliability provided by a modular 
system. A failed transmitting module would degrade system 
performance only slightly and would be easily replaceable, 

as discussed in Reference 3. 

For tracking purposes, the high altitude aircraft emits 

an omni-directional microwave beacon with a different 
frequency from the transmitting array. Interferometers on 
the ground detect phase changes in the signal received from 
the beacon. Using data from the inter- ferometers, the 
ground transmitting array can be pointed electronic-ally by 
altering the phase of each transmitting module. The system 
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Figure 10.1 Beacon Based Tracking System 


depicted in Figure 10.1 does not utilize an array, but 

rather one main trans-mitting antenna. 

The tracking systems shown above constitutes only one 
of several possible tracking system methods. Another 
possible method employs sensing devices on the aircraft. 

This system is shown in Figure 10.2 on the following page. 
These sensors detect the power received at various points on 
the aircraft. Differences in received power levels are used 
to develop signals which indicate the status of the system, 
whether the microwave beam is centered on the aircraft. A 
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simulation of this tracking system constitutes the subject 
of this Major Qualifying Project. 


10.5 Tracking System Model 

The Tracking System Model developed m this MQP is 
shown in Figure 10.3 on the following page. This Model 
consists of several major sub-systems. The first of these 
is the light source which is physically mounted to a servo 
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motor apparatus. The light and the motor simulate a 
steerable microwave beam. The Drive Circuitry block 
consists of amplifiers that provide the signals that control 
the movement of the motor. The Sensor block is made up of 
two sets of phototransistors mounted on a moveable model 
aircraft and their associated circuitry. The outputs of 
this stage 
represent the 
position of 
the light 
relative to 
the position 
of the 
sensors. 

The 

sensor 
outputs are 
sent to the 
Signal 

Processing stage where sum and difference signals are 
developed and sent to the next stage. These signals provide 
positional information regarding the airplane model with 
respect to the direction of the light beam. 

The Radio Transmitter is a 49 megahertz transmitter 
that sends the sum and difference information to the Radio 
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Figure 10.3 Model Block Diagram 
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Receiver. The Sensors, Signal Processing and Radio 
Transmitter blocks comprise the Aircraft Model. 

The signal sent by the Radio Transmitter is demodulated 
by the Radio Receiver which reconstructs the sum and 
difference signals for use in the Microprocessor and Signal 
Conversion stage. 

The Microprocessor block is comprised of an 8085A 
eight-bit microprocessor with various supporting 
peripherals. The Signal Conversion is performed by Analog 
to Digital (A/D) and Digital to Analog (D/A) converters. 
Control signals are then sent to the Drive Circuitry of the 
Motor, which completes the feedback control system. 

Upon completion of the design and construction of the 
Tracking System Model it is thoroughly tested and a controls 
analysis is completed. The testing is carried out by 
monitoring the response of the motor as the Aircraft Model 
is physically moved at different speeds and at varying 
distances from the light source. 

10.6 Summary of the Tracking System Model 

The design of the tracking system simulator proved to 
be a challenging task requiring fundamental knowledge of 
several aspects of electrical engineering. Elements of 
electronic circuit design, computer hardware design, and 
signal processing are incorporated in the Tracking System 
Model . 
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An analysis of the Tracking System Model using 
classical control theory reveals that several parameters 
exist which affect the stability and response time of the 
system. These parameters are optimized then varied to 
monitor the effects each one has on the system. Although 
the model contains several nonlinear aspects, as does every 
real system, a particularly interesting nonlinearity was 
found in the model's radio link. This nonlinearity exactly 
matches the characteristics of an ideal relay with dead 
zone, which converts a continuous input signal into a three 
level output signal. 

The operation of the Tracking System Model involves two 
modes of operation. In one of these modes, the light 
follows or tracks the motions of the model airplane. In the 
second, it performs a search pattern when the light loses 
contact with the sensors. These modes of operation are 
directly controlled by the microprocessor and can be altered 
by changing the programming of the microprocessor. The use 
of programmable memory chips makes this possible. Several 
different programs, each of which modifies the performance 
of the microprocessor in some way, are developed and 
implemented for both the track and search modes. 

Following design and testing, the operation and 
performance of the model are related to the real world 
situation. It is found that, although the Tracking System 
Model operates on a much simpler scale than an actual system 
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would, for each obstacle to achieving reliable operation 
encountered in the design of the model, there would be an 
analogous obstacle in the real system which would have to be 
dealt with accordingly. Therefore, the Tracking System 
Model provides meaningful the actual design requirements. 


»4 



11.0 Environmental Impact 


The environmental effects of operating this 
aircraft are uncertain. There is, however, evidence of 
potential harm to humans from microwave ovens which use 
the same magnetron as the microwave power station. 

Several areas of potential concern need to be studied m 
detail. The impact could take the form of possible 
heating of the surrounding air, the emission of 
microwave radiation, and the pollution or radioactive 
waste associated with the electric power generation. 

Some energy transmitted from the ground would end 
up being wasted by heating the surrounding air. There 
is an approximate 20% loss of all energy beamed from the 
transmitter (Reference 6) due to atmospheric 
attenuation. The exact long-term effects are uncertain, 
but local wildlife and human inhabitants would certainly 
be affected by changes in the climate. 

The radiation emitted from the antenna is an 
obvious source of environmental concern. Any living 
creature that might fly through the beam, whether birds 
or commercial aircraft may be exposed to dangerous 
levels of microwave radiation. There would probably be 
considerable danger to any people or animals living in 
the vicinity of the transmission source on the ground. 

Finally, there might be substantial pollution or 
radioactive waste associated with any power generator 
built for this project. While there are many existing 
power sources capable of the necessary electrical 
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output, the waste produced from the energy used would 
still be attributable to the aircraft. The lighting of a 
120 watt bulb for a year produces almost 900 pounds of 
waste products. The waste generated by this aircraft 
would dwarf this nuinbsr. 
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12.0 Cost Analysis 



The cost of building and operating this airplane 
would be tremendous. Cost of production would fall into 
two categories: building the power source and 
transmitter, and constructing the plane itself. 

The power source would be required to produce 
approximately three megawatts of electricity. Assuming 
a small nuclear reactor is the chosen method of 
providing the power, the cost would conservatively be on 
the order of several hundred million dollars. The 
transmitter for the microwave energy would have to be 
the size of a football field (assuming a wavelength of 
2.45 GHz), and would have to be able to be moved on two 
axes of rotation. Such a transmitter would cost between 
30 - 50 million dollars. 

The aircraft itself would have to be constructed 
almost entirely of composite materials which are quite 
expensive. Additionally, the reflecting surface and 
outer skin of the rectenna, and its supporting 
structure, would add substantially to the cost of 
building the plane especially is gallium-arsenide is 
used to increase efficiency (Reference 6). An estimate 
of the total airframe cost, calculated in Appendix A. 10, 
is approximately $40,237,700 (1990 dollars). 

Additional long-term operating costs would center 
on the cost of continually producing the energy required 
for the aircraft to perform its mission. 
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This study has involved the preliminary design of a 
microwave-powered, ozone-sampling aircraft for operation 
at altitudes to 100,000 feet. A first iteration of the 
design has been completed, and indicates that it may not 
presently be technically feasible to accomplish this 
mission. 

The aircraft has a gross weight of 6720 lbs. and is 
designed to carry a 1000 lb payload to an altitude of 
100,000 feet at Mach number M=0.44. The wing planform 
area of 3634 ft 2 serves as a rectifying antenna designed 
to convert microwave energy to direct current at a power 
density of 700 W/m 2 . The overall power required for the 
aircraft at the design altitude is approximately 250 
hp., while the transmitted power at the ground ranges 
from 15 to 30 MW, depending on the transmitter area. 

The aircraft utilizes both a horizontal tail and canard 
for longitudinal stability and to enhance the structural 
rigidity of the twin fuselage configuration. 

The wing structure is designed to withstand a gust- 
induced load factor n=3 at the cruise altitude, but the 
low wing loading of the aircraft makes it very sensitive 
to gusts at low altitudes. The wing structure was 
analyzed using a finite-element computational method and 
was found to be adequate at cruise altitude although 
further detailed structural design is required. Due to 
the low wing loading, the aircraft structure cannot 
withstand gusts of anticipated intensity at altitudes 



under 60,000 feet, which may induce load factors in 
excess of n=20. A gust alleviation system would have 
sensors monitoring local wind gusts and automatically 
correct for these by making adjustments to the control 
surfaces. The design of this system has not been 
conducted as part of the present project. 

The overall vehicle design has been accomplished 
using the first iteration weight and sizing estimates. 
However, a second weight iteration indicates an increase 
in gross weight, from 6720 lbs. to 15,000 lbs, and an 
associated increase in the power required to sustain 
flight. Subsequent iterations are divergent because the 
power available is a function of the wing area, which 
increases rapidly as the weight increases. This problem 
is unique to a microwave powered design. It appears 
that if the power density of 700 W/m 2 absorbed by the 
rectenna is increased, it would be possible to 
converge upon an acceptable design point. 

A beam tracking system has been developed and would 
be capable of tracking the aircraft with sufficient 
accuracy to deliver the required power. However, the 
low wing loading inherent in the design may produce 
unpredictable perturbations in the flight path which may 
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cause difficulty in tracking* 

Atmospheric absorption of the microwave radiation 
is on the order of 20% of the transmitted power, or 4 MW 
This localized heating of the air could pose a threat to 
wildlife and humans in the vicinity of the transmitter. 

A thorough environmental study should be conducted. 

This study found the cost of the airframe to be 
approximately $40 million. This value reflects tooling, 
man hours, and materials. The cost of the transmitter 
and tracking system would be 30 - 50 million dollars. 
These expenses would be excessive for the limited scope 
of the proposed mission . 
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Figure 3.3 

Airborne Transmitter Configuration 
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Figure 3.4 

Multi-Station Configuration 
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Figure 3 .5 

Single-Station Configuration 
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Figure 4.8 

First Iteration Component Weights 
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Figure 4.12 

Landing Gear Configuration 



ORIGINAL PAGE IS 
OF POOR QUALITY 


234.000 



SHEET. 













Figure 4.16 

Second Iteration Aspect Ratios and Thicknesses 
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Figure 4.17 

Second Iteration Component Weights 
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Control (MX) 
Horizontal Tal (28%) , 

Vortkal TaRt (4.7%) 


(7«e%) 


Ho ct onno (4.0%) 

tjyAwki (0.8%) 
Canard (2.4%) 


Bactranfea (24%) \ 
Landing Ooar (20%) 


Udn Wing (428%) 


v 


Motor* (7.8%) 
•artei . 

PTOpWirl 




Tqylood (8.7%) 
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Figure 5,1.2 Sample Pressure Distribution 
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Figure 5.1.4 L, P, Q Target Pressure Distributions 
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Figure 5.1.5 P-Series Pressure Distribution (inviscid analysis) 


Airfoil PI 
design point 
cl c 0.7238 
alpha-0. 314° 







Figure 5.1.7 Q-Series Pressure Distribution (inviscid analysis) 


Airfoil Q1 
design point 
Cj-0.740 




Target Pressure Distributions 








Figure 5.1.11 L3 Inviscid Pressure Distribution, Geometry 


Airfoil L3 
design point 
c =0.399 
alpha=0. 35° 




iscid analysis) 

Attack 



Section cl ve 
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Figure 5.1.14 

CP1NVRS.F0R 
Trailing Edge Closure 


Extrapolated 
Blunt Trailing Edge 









Airfoil Coordinates 

Zero Thickness Trailing Edge 
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Analysis 




Figure 5.1.19 LI, NACA 4415 Drag Comparison 

cd versus cl 









Figure 5.1.22 LI, NACA 4415 Lift Comparison 

Section d versus Angle of Attack 



alpha (dag.) 

D LI + 2UCA 


Figure 5.1.23 L2, NACA4412 Lift Comparison 

Section cl versus Angle of Attack 



alpha (deg.) 

+ KJLCA 4412 


D 


L2 


FiEure 5.1.24 L3, NACA 4409 Lift Comparison 

cl versus alpha 

flnmoirlaon Of L3 With 1UQU409 



alpha 

□ u + HACA. 4409 


cl vs % span 

Are. lift coefficient of 


Figure 5.2.1 



span loacation 


Figure 5.2.2 

Airfoil Lift Curve Equations 


Airfoil LI: 

c 1 = 7.02 a + .673 

Airfoil L2: 

c 1 = 6.80 a + .504 

Airfoil L3 : 

c 1 = 6.70 a + .359 


Figure 5.2.3 

Variation of Lift Curve Slope 



lift 

curve slope: 

0<b< . 6 

m = 

-.3667b + 7.02 

. 6<b< . 9 

m = 

-.3333b + 6.99 

.9<b< 1 

m = 

constant = 6.69 per radian 


Figure 5.2.4 

Variation of Zero Lift Angle of Attack 


zero 

angle lift coefficient: 

0<b< . 6 

c. = -.2813b + .673 

1 Q = 0 

.6<b<.9 

c. = -.4852b + .7953 

la-0 

. 9<b< 1 

c. = constant = .3586 

1 O K 0 








Twist Distribution 


Figure 5- 2 
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TWIST 0.0 DEG. 
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Figure 5.2.8 
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Figure 5.2.13 

Untwisted Plan-form 
C L = 0.4237 
alpha = -1.0° 




Figure 5.2.14 


Untwisted Planform 
C = 0.600 
alpha = 0.375° 
Design Point 



Figure 5.2 

Untwi sted 
C, = 0.814 
alpha = 2. 
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Untwisted Planform 
C, = 1.2 
alpha = 5.0° 
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Figure 5.2.18 

C, = 0.6 
alpha = 0.2° 
Twisted Planform 
Design Point 










FIGURE 6.1.1 
Spanwise & Chordwise 
Pressure Distribution 
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Figure 6.1.4 
Model #1 
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Figure 6.1.5 
Model #1 
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Figure 6.1.7 
Model #1 


(V cn 

< CO 

m oi 

• ** 

w 

CD 

> 

0) c 

Z ID 

< Q 


CVJ 


ID 

> 

< 


M/) 1 I 
^ UJ CL CT 
..CCIDUJX 
CD H H H X 
n [/) l/] H U) 


O 

*•« 

O (D 
in ^ 'H 
ID U I CVJ 
_) 8 8 
D X 

Q ID Z X 
w _j 2 2 

x id cn in 







ID 

CVJ 

N 













D 

in 





in 









D 

CVJ 

cn 


ID 

N 

I s - 









M 

• 

• 

• 

•<T 

cn 

CVJ 






QD 


l 

I 

OJ 

I s - 

CVJ 


• 

• 

s 

CVJ 

N 

o 



o 



cn 

I s - 

nj 

CVJ 

N 

CVJ 

01 

m 

o 

o 

o 

o 

U) 


UJ 

CD 



in 

CD 


CVJ 

(0 

o 

^-4 CD 

a! 

t-4 

ID 

cn 

in 

1 

i 

I 

CVJ 

in 

U) 

▼4 

-*-1 

CVJ 

1 8 

n 

1 

B 

5 

M 

1 

i 

0 

K 

8 

s 

c 

8 

8 

H 




O 

u 










w 




H 

ID 










XV 

DI 

XF 

U- 

X 

ZF 

c r 

cr 

a 

< 

ID 

u 

n 

ID 

tL 

ID 

I 

w 


I 


I 

i 


l 

1 


I 


i 

I 

i 

1 

I 

i 

t 

! 

J 


! 



ORHSNL*,!. PAGE J3 
Or PCG8 QUALITY 


j 

i 


i 

I 

1 

I 

J 


cu 

ID 



L 


Left I beam Weight = 2050 lbs. BAR 


Figure 6.1.8 
Model #1 
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Figure 6.1.9 
Model #1 
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Figure 6.1.10 
Model #1 
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Figure 6.1.11 
Model #1 
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Figure 6.1.12 
Model #1 






Figure 6.1.13 
I Beam Properties 


Section Thick 
Front 
I Beam 
1 0.1 

2 0.09 

3 0.08 

4 0.07 

5 0.06 

6 0.063 

7 0.046 

8 0.0426 

9 0.04 

10 0.04 

20 0.04 


Thick 

Thick 

Middle 

Rear 

I Beam 

I Beam 

0.11 

0.105 

0.098 

0.094 

0.087 

0.083 

0.079 

0.075 

0.07 

0.065 

0.058 

0.055 

0.047 

0.046 

0.043 

0.042 

0.04 

0.04 

0.04 

0.04 

0.04 

0.04 



Figure 6.1.14 
TABLE 6.1.6 


SECTION 

AREA THICK 

BASE 

Ixx 

iyy 


FRONT FRONT 

FRONT 




I BEAM I BEAM 

I BEAM 



1 

2.16 0.5 

4.32 

0.045 

3.36 

2 

1.94 0.45 

4.32 

0.0328 

3.02 

3 

1.728 0.4 

4.32 

0.023 

2.69 

4 

1.52 0.35 

4.32 

0.015 

2.35 

5 

1.32 0.3 

4.32 

0.00972 

2.02 

6 

1.14 0.26 

4.32 

0.00633 

1.75 

7 

0.972 0.225 

4.32 

0.0041 

1.51 

8 

0.918 0.21 

4.32 

0.00333 

1.41 

9 

0.864 0.2 

4.32 

0.00288 

1.34 

10 

0.864 0.2 

4.32 

0.00288 

1.34 

20 

0.864 0.2 

4.32 

0.00288 

1.34 


Figure 6.1.15 





TABLE 6.1. 

.7 



SECTION 

AREA THICK 

BASE 

Ixx 

iyy 


MIDDLE MIDDLE 

MIDDLE 




I BEAM I BEAM 

I BEAM 



1 

1.872 0.5 

3.74 

0.039 

2.18 

2 

1.68 0.45 

3.74 

0.028 

1.96 

3 

1.49 0.4 

3.74 

0.02 

1.74 

4 

1.327 0.355 

3.74 

0.014 

1.55 

5 

1.164 0.31 

3.74 

0.00928 

1.35 

6 

0.995 0.27 

3.74 

0.00613 

1.18 

7 

0.825 0.22 

3.74 

0.00332 

0.96 

8 

0.7725 0.207 

3.74 

0.00276 

0.9 

9 

0.72 0.19 

3.74 

0.00214 

0.83 

10 

0.72 0.19 

3.74 

0.00214 

0.83 

20 

0.72 0.19 

3.74 

0.00214 

0.83 


Figure 6.1.16 





TABLE 6.1 

.8 



SECTION 

AREA THICK 

BASE 

Ixx 

Iyy 


REAR REAR 

REAR 




I BEAM I BEAM 

I BEAM 



1 

1.584 0.5 

3.17 

0.033 

1.33 

2 

1.418 0.45 

3.17 

0.024 

1.195 

3 

1.253 0.4 

3.17 

0.017 

1.06 

4 

1.131 0.357 

3.17 

0.012 

0.948 

5 

1.01 0.319 

3.17 

0.00858 

0.847 

6 

0.844 0.266 

3.17 

0.00497 

0.706 

7 

0.677 0.214 

3.17 

0.00259 

0.568 

8 

0.626 0.197 

3.17 

0.00202 

0.523 

9 

0.576 0.182 

3.17 

0.00159 

0.483 

10 

0.576 0.182 

3.17 

0.00159 

0.483 

20 

0.576 0.182 

3.17 

0.00159 

0.483 


Figure 6.1.17 
Model n 
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Figure 6.1.18 
Model H2 
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Figure 6 . 1 .41 
Model #3 
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Figure 6.1.42 
Model #3 
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Figure 6.2.1 


The frequency value in each mode for 
preliminary wing structural design is shown below: 

EIGENVALUE (NATURAL FREQUENCY) 

MODE FREQUENCY (CYCLE/SEC) 

1 0.858065 

2 4.048 

3 4.843 

4 11.931 

5 13.368 



FIGURE 6.2.2 

Mode 1 Preliminary Design 
1st Bending 
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Mode 1 Preliminary Design 
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FIGURE 6.2.4 

Mode 1 Preliminary Design 
1st Bending 
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FIGURE 6.2.5 

Mode 2 Preliminary Design 
1st Torsion 
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FIGURE 6t2.6 

Mode 2 Preliminary Design 
1st Torsion 
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FIGURE 6.2. 7 

Mode 2 Preliminary Design 
1st Torsion 
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FIGURE 6.2.8 

Mode 3 Preliminary Design 
2nd Bending 
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FIGURE 6.2.9 

Mode 3 Preliminary Design 
2nd Bending 
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FIGURE 6.2.10 

Mode 3 Preliminary Design 

2nd Bending 
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FIGURE 6.2.11 

Mode 4 Preliminary Design 
Combined mode 
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FIGURE 6.2.1? 

Mode 4 Preliminary Design 
Combined mode 
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FIGURE 6.2.13 

Mode 4 Preliminary Design 
Combined mode 
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FIGURE 6.2.14 

Mode 5 Preliminary Design 

Combined mode 
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FIGURE 6.2.15 

Mode 5 Preliminary Design 
Combined mode 
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FIGURE 6.2.16 

Mode 5 Preliminary Design 
Combined mode 
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Figure 
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The frequency value 


in each mode 


for the revised 


wing : 


EIGENVALUE (NATURAL FREQUENCY) 

MODE FREQUENCY (CYCLE/SEC) 

x 0.937703 

2 5.27937 

3 7.20949 

4 14.0685 

5 21.3167 

6 26.8827 

7 35.9328 

8 43.5590 

9 52.0220 

in 63 . 4687 
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FIGURE 6.2.18 

Mode 1 Revised Design 

1st Bending 
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FIGURE 6.2.19 

Mode 1 Revised Design 

1st Bending 
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FIGURE _6 . 2 . 20 
Mode 1 Revised Design 
1st Bending 
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FIGURE 6.2.21 

Mode 2 Revised Design 

2nd Bending 
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FIGURE 6.2.22 

Mode 2 Revised Design 

2nd Bending 


CXI o 
< cn 
<n 01 


OD 

to w 

> 

tO 31 
2 UJ 
< U. 


*H 0) 


01 O) 
N ® 

OJ tu 

• cu in ** 

j i c n i 

niro 
to uj u m x 

H h H I 2 

QWHUQ 


0) 

in _ 

*4 N 01 

i : 

H 
CO 

> H U. IL 
X D x > 


o o 


I fl I 


0 

01 
ID 

I 


u. 

N 



AMSYS 



10 . 


FIGURE 6.2.23 

Mode 2 Revised Design 

2nd Bending 





FIGURE 6.2.24 

Mode 3 Revised Design 

1st Torsion 
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FIGURE 6.2.25 

Mode 3 Revised Design 

1st Torsion 
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FIGURE 6.2.26 

Mode 3 Revised Design 

1st Torsion 
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FIGURE 6.2.27 
Mode 4 Revised Design 
3rd Bending 
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FIGURE 6.2.28 

Mode 4 Revised Design 
3rd Bending 
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FIGURE 6.2.29 

Mode 4 Revised Design 

3rd Bending 



CO 






10 

N 





o 



O 




CU 

cn 

• 

o 




in o 

o 

01 

Ti 

T -H 


N 01 


10 

1 

: ii 

B 

1 D S 

I 

c 

a 

tr a 


h* 



UJ 

\u LU 

■/ 

w 



H 

H 2C 


> H lL 

lL 

u. 

in 

►s Ll 

n 

>- Q X 

>■ 

N 



1W1NG3 



FIGURE 6.2.31 

Mode 5 Revised Design 

2nd Torsion 
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FIGURE 6.2.32 
Mode 5 Revised Design 
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The frequency value in each 

mode for the revised 


wing : 




EIGENVALUE (NATURAL FREQUENCY) 


MODE 

FREQUENCY (CYCLE/SEC) 


1 

1.54863 


2 

5.88586 


3 

7.17381 


4 

12.6177 

— 

5 

20 . 9103 


6 

25.4295 


7 

35 . 9099 


8 

37 . 5427 


9 

47.5733 


10 

56 . 8922 





FIGURE 6.2.34 

Mode 1 Revised Design w/Brace 
1st Bending 
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FIGURE 6.2.35 

Mode 1 Revised Design w/Brace 
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FIGURE 6.2.36 

Mode 1 Revised Design w/Brace 
1st Bending 
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FIGURE 6.2.37 

Mode 2 Revised Design w/Brace 
2nd Bending 
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FIGURE 6.2.38 

Mode 2 Revised Design w/Brace 
2nd Bending 
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Mode 2 Revised Design w/Brace 
2nd Bending 
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FIGURE 6.2.40 
Mode 3 Revised Des 
1st Torsion 
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FIGURE 6.2.41 

Mode 3 Revised Design w/Brace 
1st Torsion 
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FIGURE 6.2.42 

Mode 3 Revised Design w/Brace 
1st Torsion 
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FIGURE 6.2.43 

Mode 4 Revised Design w/Brace 
3rd Bending 
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FIGURE 6.2.44 

Mode 4 Revised Design w/Brace 
3rd Bending 
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Mode 4 Revised Design w/Brace 
3rd Bending 
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FIGURE 6.2.47 

Mode 5 Revised Design w/Brace 
2nd Torsion 
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FIGURE 6.2.48 

Mode 5 Revised Design w/Braee 
2nd Torsion 
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Static Stability Analysis Data 

• nescrintion : 

Value: 

o y mu'w' ±_ 
C 

wing mean aerodynamic chord 

15.36 ft. 

S w 

wing planform area 

3594.3 sq. 


wing aspect ratio 

15.23 

s t 

tail planform area 

520 sq. ft 

AR t 

tail aspect ratio 

5.2 

S c 

canard planform area 

470 sq. ft 

^c 

canard aspect ratio 

4.7 

L t 

CG to tail AC distance 

36.50 ft. 

L c 

CG to canard AC distance 

23.66 ft. 

y 

CG to wing AC distance 

3.84 ft. 

ac 

y 

CG to neutral point dist. 

3.05 ft. 

eg 

n 

wing/tail effective vel. ratio 

1 

c L o» 

wing lift coeff. at 0 AOA 

0.302 

Cm ac" 

wing moment coeff . at 0 AOA 

-0.0596 


deader wing section lift curve slope 


27T 
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Figure 7.2.3 

Short Period Mode ■ 100,000 ft 
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Figure 7.2.4 

Long Period Mode - 100,000 ft 
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Figure 7.2.5 

Combined Mode - 100,000 ft 
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Figure 7.2.6 

Short Period Mode - 25,000 ft 




Time (seconds) 



it, (!l-axis) vs 1NDEP.VAR (X-axis) 


Figure 7.2.7 

Long Period Mode - 25,000 ft 
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Figure 7.2.8 

Combined Mode - 25,000 ft 
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Figure 8.2 

Component Locations and Mass Properties 
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Power Required 

Wings Level Flight 


Figure 9.2.2 
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Power Available and Required 

Altitudo - Sea Level (level flight) 



Mach Number 

Power Required Power Available 



Power Available and Required 

Altitude - 25,000 foot (level flight) 



Mach Number 

Power Required Power Available 



Power Available and Required 

Altitude - 50,000 foot (level flight) 



Mach Number 

Power Required * Power Available 



Power Available and Required 

Altitude - 75,000 foot (level flight) 



Power Required Power Available 



OOG 


O 

o 


o 

o 

CO 


o 

o 

CM 


o 

o 


(dq) aoAVOcJ 






Figure 9.2.8 

Angle of Incidence 


e 

bank angle 

<P 

beam angle 

6 

incidence 

\ \ 

angle 


altitude 

i r 

flight 

radius 

\ 1 / 


w 




6 = <p + © 


e = tan-l (R/h) 


<p = tan-1 ( V2 /(g 

* R)) 

tan-l (R/h) + tan-1 (V2 /(g * R) ) 




Angle of Incidence 


Figure 9.2.9 


Flight Radius Increases 


Lim © = 90* 
R- oo 


Lim tp = 0* 

R- oo 


lim B = 90* 
R- oo 
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Figure 9.2.10 


Flight Radius Decreases 


lim © = 


0* 


v) l 







Power Required 



(dti) asAiOd 




Net Power Available 

Altitude — SO, 000 feet 


Figure 9.2.13 



(dn) 


Mach Numbor 

Bank - 10 Bank = 15 Bank 


Power Required 



Mach Number 





Net Power Available 

Altitude - 75,000 feet 
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n 



(dL-ci) aQM OcJ 


Mach Number 

Bank « 10 Bank - 15 Bank « 20 



Power Required 



Mach Number 




Net. Power Available 

Altitude - 100,000 feet 


Figure 9.2.17 
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Figure 9.2.19 


Transmitted Power Distribution 



Pd power 

density 

At transmitter 
area 

Dt transmitter 
diameter 
Dr receiver 
diameter 
BVI\ \ beam 
h altitude 


BW*j = 50 * wavelength * (27T/360V 
Dt 

Dt = 50 * wavelength * (27r/360^ 
BW^ 

Geometric Relationship 


tan ( BW^/2 ) = ^Dr/h 
BW^ = 2 * tan-1 ( ^ Dr /h) 


Solving for Dt in terms of Dr : 


Dt = ?5 * wavelength * f27T/360l « 10 * wavelength * 7T 
tan-1 ( IjDr/h) 36 * © 


Figure 9.2.20 


Parabolic Power Density Distribution 

Y fraction of Pd 
X ratio of 
span to Dr 

Pa 



Parabolic Distribution Described by : 

r = —3$X2 h 


Figure 9.2.21 


Total Transmitted Power Calculation 

(at cruising altitude) © BW^ 


pt = Pd * wavelength * h2 (9.9) 
At * n 


tan (6/2) ©(rad) Dt(m) At(m2) Pd(W/m2) P(MW) 


1.25E-3 2.5E-3 42.761 1436.1 123.74 15.016 


1.50e-3 3.0E-3 35.624 997.27 1002.3 17.514 


1.75E-3 3.5E-3 30.540 732.69 899.27 21.388 


2.00E-3 4.0E-3 26.725 560.97 843.02 26.188 


2.50E-3 5.0E-3 21.380 359.01 785.26 38.117 







Rate of Climb vs Altitude 

Wings Level Flight 


Figure 9.3.1 
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Rate of Climb vs Mach Number 


Figure 9.3.2 
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Rate of Climb vs Mach Number 

Altitude - 90.000 feet 


Figure 9.3.4 
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Rate of Climb vs Mach Number 


Figure 9.3.5 
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Rate of Climb vs Mach Number 

Altitude - 100,000 feet 


Figure 9.3.6 
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Figure 9.3.7 


Time to Climb 
Wings Level Flight 


Time 

altitude dh/dt(f/s) dh/dt(f/s)-l seconds minutes 


Sea Level 

to 


25,000 

ft 

25,000 

ft 

to 


50,000 

ft 

50,000 

ft 

to 


75,000 

ft 

75,000 

ft 

to 


100,000 

ft 


12.381 

0.080772 

2019.3 

33.655 

11.571 

0.086424 

2160.6 

36.010 

9.5839 

0.10434 

2608.5 

43.475 

5.1744 

0.19315 

4828.7 

80.478 
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Bank Angle vs Altitude 


Figure 9.4.5 
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Flight Radius vs Altitude 


Figure 9.4.6 
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Climb Rate vs Altitude 


Figure 9.4.7 
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Altitude vs Flight Radius 


Figure 9.4.8 
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EFFECT OF WIND GUST 

ALTITUDE - SEA LEVEL 


Figure 9.5.1 


V-n Diagram 
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EFFECT OF WIND GUST 

ALTITUDE - 25,000 FEET 


Figure 9.5.2 
V-n Diagram 
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EFFECT OF WIND GUST 
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EFFECT OF WIND GUST 

ALTITUDE - 50,000 FEET 
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EFFECT OF WIND GUST 

ALTITUDE - 100,000 FEET 


Figure 9.5.10 
V-n Diagram 
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Figure 9.5.11 
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Figure 9.5.13 
V-n Diagram 
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